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^Jl'his article considers recent advances in the investigation of the thermal and magnetic properties of 
^integrable spin ladder models and their applicability to the physics of strong coupling ladder com- 
QJpounds. For this class of compounds the rung coupling Jj_ is much stronger than the coupling J« 
' yllpng the ladder legs. The ground state properties of the integrable two-leg spin-i and the mixed 
^y^spin-(^, 1) ladder models at zero temperature are analysed by means of the Thermodynamic Bethe 
K^A'nsatz (TBA). Solving the TBA equations yields exact results for the critical fields and critical 
'^behaviour. The thermal and magnetic properties of the models are discussed in terms of the recently 
introduced High Temperature Expansion (HTE) method, which is reviewed in detail. In the strong 
Q?c oupling region the integrable spin-^ ladder model exhibits three quantum phases: (i) a gapped phase 
^■J iri the regime H < H c \ = J± — 4 Jii , (ii) a fully polarized phase for H > H c 2 = J± + 4 Ju , and (iii) a 
vNLuttinger liquid magnetic phase in the regime H c i < H < H c 2- The critical behaviour in the vicinity 
H_f the critical points H c i and H c 2 is of Pokrovsky-Talapov type. The temperature-dependent thermal 
l/~)and magnetic properties are directly evaluated from the exact free energy expression and compared 
C _ X o known experimental results for the strong coupling ladder compounds (5IAP)2CuBr4 ■ 2H2O, 
■^JC»2(C5Hi2N2)2Cl4, (C 5 Hi2N) 2 CuBr4, BIP-BNO and [Cu2(C20 2 )(CioH 8 N 2 )2)](N03)2. Similar 
jH^nalysis of the mixed spin-(^,l) ladder model reveals a rich phase diagram, with a -| and a full 

Saturation magnetization plateau within the strong antiferromagnetic rung coupling regime. For 
/eak rung coupling, the fractional magnetization plateau is diminished and a new quantum phase 
I ' 

Or 



Ot 

ot 

•l>: 
X 



Murray.Batchelor@anu.edu.au 
xwcl05@rsphysse.anu.edu.au 
oelkers@maths.uq.edu.au 
zengo_tsuboi@prcf.okayama.jp 



Advances in Physics ISSN 0001-8732 print/ ISSN 1460-6976 online ©2005 Taylor & Francis Ltd 
http: //www. tandf.co.uk/journals 
DOI: 10. 1080/00018730xxxxxxxxxxxx 



2 



M.T. Batchelor, X.-W. Guan, N. Oelkers and Z. Tsuboi 



transition occurs. The phase diagram can be directly deduced from the magnetization curve obtained 
from the exact result derived from the TBA and HTE. The results are applied to the mixed ferri- 
magnetic ladder compound PNNBNO. The thermodynamics of the spin-orbital model with different 
single-ion anisotropics is also investigated. For this model single-ion anisotropy can trigger differ- 
ent quantum phase transitions within the spin and orbital degrees of freedom, with magnetization 
plateaux arising from different spin and orbit Lande g-factors. 



1 Introduction 

The subject of integrable or exactly solved models in statistical mechanics 
has inspired a number of profound developments in both physics and mathe- 
matics. There are several key solvable models in low dimensions, the analysis 
of which has developed to an extraordinary level of mathematical sophistica- 
tion. This article is concerned with the calculation of physical properties of 
the integrable quantum spin ladders and their applicability to experimental 
ladder compounds. By integrable it is meant that the ladder models are ex- 
actly solvable in the Yang-Baxter sense [IHE] > with infinitely many conserved 
quantities and an underlying Bethe Ansatz solution. Most importantly, this 
entails access to the thermodynamic properties [4, 5J. In this way the theory 
of integrable models, which has developed for over 75 years since Bethe's [6l[7] 
pioneering solution of the su(2) Heisenberg chain, has recently been applied [8] 
to make direct contact with experiments on spin ladder compounds. The in- 
tegrable ladder models, for example, the two-leg spin-^ ladder and the mixed 
spin-(^, 1) ladder, are variants of the su(A) and su(6) models based on the 
general family of integrable su(N) permutator models fl^ [2|, fTO], 1 1 1 j . Their 
two-dimensional classical lattice model counterparts are the iV-state An_\ 
vertex models [EJ[E3], which are special cases of the Perk-Schultz model [13] . 

1.1 Physical motivation 

Over the last few decades there has been considerable theoretical and ex- 
perimental interest in low-dimensional quantum spin systems. One focus of 
attention is the spin liquid behaviour. It is well established that the critical 
behaviour of the one-dimensional Heisenberg spin chain in the absence of frus- 
tration depends on the spin magnitude S, i.e., on whether 2S is an even or 
odd integer. In one dimension the quantum fluctuations are so strong that 
the long range Neel order no longer survives. For 25 odd the spin chain has 
low-lying gapless excitations with a power law spin-spin correlation [15]. On 
the other hand, for 2S even the spin-spin correlations decay exponentially due 
to the Haldane gap [161117]. In two-dimensional quantum spin systems the 
quantum fluctuations compete with the thermal fluctuations for domination 
of the critical behaviour at low temperatures. Here long range Neel order is 
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possibly stable against quantum fluctuations. 

Spin ladder systems, which intermediate between one-dimensional and two- 
dimensional spin systems, exhibit rich and novel phase transitions and a di- 
verse range of quantum magnetic effects at low temperatures jl8 t ll9 1 [2"0" l l21j. 
The long range Neel order phase emerges when the interaction between the 
adjacent two-leg ladders is greater than some critical value [22]. However, if the 
two-leg ladder structures are well isolated from each other, the intra- and inter- 
chain interactions induce rich spin liquid phases, such as gapped dimerised and 
gapless phases. The existence of a spin gap, magnetization plateaux, quantum 
critical points and superconductivity under hole doping are examples of key 
physical properties observed in the ladder compounds. The physics of the 
spin ladder systems reveals a significant pairing effect, namely the Heisen- 
berg spin ladders have an energy gap in the spin excitation spectrum if the 
number of ladder legs is even, whereas there is no gap for an odd number of 
legs |18 [ ll9 [ l20]. This similarity with the above mentioned Haldane conjecture 
for spin-S chains arises from the fact that the ladder model can be mapped 
onto a one-dimensional counterpart with on-site anisotropy or next-nearest- 
neighbour interaction, at least in the strong coupling limit [23U24j . In addition, 
the ladder systems have been argued to provide simple models relevant to the 
mechanism for high temperature superconductivity in layered cuprate super- 
conductors [20]. The challenge for theorists is to calculate and to predict the 
properties of the ladder systems. 

1.2 Experiments and theories for ladders 

The experimental study of ladder systems was initiated from the observation 
that the compound (VO)2P207 has a singlet ground state with an energy gap 
in the spectrum of spin excitations [25j . Subsequently, the theoretical study 
of ladder systems began when the hole-doped two-leg t—J ladder was found 
to exhibit a finite spin gap [26]. In this model, at each site the ions are either 
in a spin-^ Cu 2+ state or a Cu 3+ hole state. Hole ions are bound along rungs 
so that the spin-^ Cu 2+ ions form rung singlets. This configuration leads to 
superconducting pairing correlations [26|,l20j . In this context, there have been 
many theoretical papers focussing on hole-doped ladders, see Refs. |20[I27]. 
Experimental superconducting ladder compounds, e.g., (La,Sr,Ca)i4Cu2404i 
[28|l29] and other doped ladder compounds, e.g., LaCu02.5 with Sr doping [30j 
127] . provide experimental evidence of superconductivity. In Cu-based ladders, 
the hole doping provides a "d-wave" channel for a superconducting phase 
|26U20j . In this article we focus on undoped or pure spin ladder materials for 
comparisons with the theoretical results. 

The Heisenberg antiferromagnetic spin-^ ladder model was first proposed 
in Ref. [18]. It is known that the Heisenberg coupling J± along the rungs can 
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open a gapped phase for arbitrary coupling strength. If J± = 0, the ladder 
decouples into two non-interacting spin-^ Heisenberg chains, with no gap to 
spin excitations [15]. The decoupled chains are in critical states exhibiting 
no long range order. In the presence of weak interchain interactions (small 
J_i_), the lowest spin excitation can be separated from the ground state by an 
energy gap of magnitude A « Jj_/2 • Essentially, the weakly 

coupled ladder models may exhibit antiferromagnetic long-range order at finite 
temperatures. For example, the trellis layer structure of the Sr(Cui_ x Zn x .)203 
lattice with nonmagnetic isoelectronic Zn +2 doping |34U35| was found to induce 
long-range Neel order. It was verified experimentally that the undoped two- 
leg ladder compound SrCu203 exhibits a spin gap, but the three-leg ladder 
compound S^CusOs does not [36][371[38]. In the strong coupling limit, the 
ground state of the two- leg ladder compound is dimerised along the rungs. The 
intrachain interaction between neighbouring magnetic moments is quenched 
by the formation of a rung singlet state. In order to create a spin excitation 
the rung singlet must be replaced by a rung triplet. Breaking the rung singlet 
causes an energy gap A J± — Jji which separates the lowest excited state 
from the dimerised ground state [20 . 22. 39] . Here Jy is the (intrachain) coupling 
along the ladder legs. The number of legs significantly affects the magnitude of 
the spin gap. For an even number of spin-^ ladder legs the gap decreases as the 
number of legs grows |36|I37|. as one would expect as the ladder becomes more 
like a two-dimensional antiferromagnet. The ground state is a spin singlet 
with zero total spin. For an odd number of legs with strong rung coupling, 
the ground state is a magnetic degenerate state with total spin S = i, with 
algebraically decaying spin-spin correlations and no gap to spin excitations. 
The qualitatively different behaviour between even- and odd-leg spin ladders 
has been discussed at length plH0lHTl[20] . 

Of particular interest are the strongly coupled ladders as their energy 
gaps are experimentally accessible. Many organic spin ladders with strong 
rung coupling have been synthesised. In particular, (5IAP)2CuBr4-2H20 [32] 
Cu 2 (C 5 Hi2N2)2Cl4 PIOTSES] and (C 5 Hi2N) 2 CuBr4 [47] have a well isolated 
ladder structure. The ladder structure for the compound Cu2(C5Hi2N2)2Ci4 
is shown in figure [TJ This compound is made up of effective two- leg spin-^ lad- 
ders with strong antiferromagnetic rung interaction through superexchange 
between adjacent molecular units. The interchain coupling between Cu-Cu 
ions is parametrized by the rung coupling J±. The dimeric units along each 
rung are weakly linked to the nearest neighbour ones through hydrogen bonds 
along the [101] chain direction. The exchange between two neighbouring rungs 
is much weaker than the intradimer path. This intrachain coupling is denoted 
by the parameter Jm. It is believed that the family of compounds KCuCbj, 
TlCuCi3, NH4CUCI3 and KCuBr3 also exhibit such double-spin-chain struc- 
ture PlHgi^l^l^ESllHlEaiSlEZlESlEgiEDlEIlEZ]. These double-chain 
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Figure 1. The schematic structure of the two-leg spin ladder compound Cu2(C5Hi2N2)2Cl4. The 
Cu-Cu binuclcar units along the chain direction [101] form the intrachain interaction. The 
supcrcxchangc paths between two Cu ions along each rung form the interchain interaction. 

Reproduced from Ref. |46| . 

compounds have been studied extensively from both the theoretical and ex- 
perimental points of view. 

The Heisenberg ladder model, with antiferromagnetic spin-spin interactions 
along the legs and rungs, describes these compounds rather well. The strong 
coupling compounds have a singlet ground state for magnetic fields less than 
a critical value H c \. The theoretical prediction for the energy gap is A ~ 
J_i_ — Jii 163 ^ 164 ^ 120]. which is in excellent agreement with the experimental 
results, obtained from Nuclear Magnetic Resonance measurements [44], along 
with high field magnetic properties [59. 42. 65 . 50]. In addition, the susceptibility 
X oc e~ A//T / y/T with energy gap A obtained from the quantum transfer matrix 
method [66] generally provides a good fit for the strongly coupled two-leg 
ladders. For magnetic field H > H c \ the triplet rung state enters into the 
ground state, with a gapless Luttinger phase in the regime H c \ < H < H C 2. 
The critical field H C 2 indicates a quantum phase transition separating the 
Luttinger liquid magnetic phase from a fully polarized state. It was suggested 
[67] that the gapless phase can be viewed as a Bose-Einstein condensation of 
magnons [68I69U70II71P72P73] . The quantum fluctuation of low- lying excitations 
for these quantum spin systems are characterized by spin waves which reveal 
bosonic signature. In this way the strongly coupled spin ladder with magnon 
excitations in high magnetic field can be mapped onto a one-dimensional XX Z 
Heisenberg model with an effective magnetic field. From this mapping, the 
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critical fields are given by H c \ « Jj_ — J« and H C 2 ~ J± + 2Jy [E2[M|> which 
agree rather well with experiments [20] . 

Considerable interest has been generated by fractional magnetization 
plateaux, which have been found in Shastry-Sutherland dimer systems [7U 
[751 [76] an d some spin organic ferrimagnetic ladders [77J [7HJ [79] . Theoreti- 
cal studies and numerical results suggest that magnetization plateaux ex- 
ist in a mixed spin ladder [80l[8TJ[82] and the mixed spin-(^,l) chains 
pl[8il[85l[86l[871[HH|[89l[90] . It was concluded that for certain strong rung cou- 
pling values, there exist magnetization plateaux at M z = 0.5 and at M z = 1, 
but no plateaux for negative (ferromagnetic) rung coupling values. In the 
strong coupling limit, the mixed spin ladder can be viewed as a ferrimagnetic 
chain which displays unusual quantum effects at low temperatures. 

1.3 Integrable ladder models 

Haldane's conjectured universal properties of the spin-S" Heisenberg antiferro- 
magnet depending on whether 2S is odd or even have been recognized in many 
disordered spin liquid systems, such as in spin ladders, mixed spin ladders and 
alternating spin chains. The spin liquid systems have a finite correlation length 
and a spin gap at zero temperature. The ground state is a spin singlet. The 
lowest spin triplet excitations have an energy gap due to the finite correlation 
length. In this context, Nersesyan and Tsvelik [91] predicted a different gapped 
phase for a generalised two-leg spin-Jj ladder model with a strong four-spin in- 
teraction term (the biquadratic interaction strength \ V\ is very large [91]). This 
model includes the conventional two- leg Heisenberg spin-^ ladder and the in- 
tegrable spin ladder model based on su(4) algebra symmetry [S2] as special 
cases, i.e., V = and |V| = 4Jy, respectively. In the case of weak biquadratic 
interaction (|V| very small), the generalised Heisenberg ladder model has the 
same Haldane spin liquid phase as that of the conventional spin ladder. How- 
ever, for strong enough biquadratic interaction this generalized ladder model 
exhibits a spontaneous dimerization phase which is different from the Haldane 
spin liquid gapped phase. The physical importance of the model has been 
subsequently noted [93J. 

From the theoretical perspective, these ladder models cannot be solved ex- 
actly, in contrast to the spin-^ Heisenberg chain. Thus their thermal and 
magnetic properties are not accessible via exact methods. In order to get a 
handle on the physical properties, some integrable extensions of the Heisen- 
berg ladders have been constructed [MllMll^lMllWllMllMl flT^fll^flTOflWl 
^l^^^j^jj^^n]^E^^SI^I^IW- These models include an 
integrable ladder based on the su(4) algebra [92] and integrable ra-leg ladder 
models [96,97,98,99,114,115] among a number of variants. At first glance, due 
to the complicated spin-spin interactions along the legs and rungs, it appears 
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that none of these models can be used to predict physical properties for direct 
comparison with experimental data. Indeed, the application of the integrable 
ladder models to ladder materials should be taken with caution. In this article 
we focus on the calculation and study of the exact thermal and magnetic prop- 
erties of the integrable ladder models and attempt to fit these properties to 
the experimental data for some known strong coupling ladder compounds by 
setting aside the issues regarding the applicability of the integrable Hamilto- 
nian to those materials, in particular the effect of the biquadratic interaction 
term. We demonstrate that under certain conditions, namely 

Ji » J|| and K B T > J\\ (1) 

the integrable ladder models may be used to study the thermal and magnetic 
properties of ladder materials with strong rung coupling. However, at low tem- 
peratures the biquadratic interaction term in the integrable models may lead 
to discrepancies with the physics of the ladder materials. These discrepancies 
are minimal if J± 3> Jji. We point out that the good agreement between the 
theoretical predictions and the experimental results for macroscopic quantities 
such as the magnetization, the susceptibility and the specific heat do not im- 
ply that the materials considered have the integrable Hamiltonian structure. 
For determining the precise structure of the materials more physical analysis, 
such as dynamical structure and dispersion in low-lying excitations, should 
be carried out in addition to the study of the thermal and magnetic proper- 
ties. At high temperatures, the thermodynamics of the compounds can also 
be modelled by a rung-dimer system [185] . However, the dimer model gives 
poor agreement with experimental results at low temperatures. For the dimer 
model, it appears necessary to include Dzyaloshinskii-Moriya interaction terms 
in order to better explain the experimental behaviour [186] , 

The integrable ladder model based on the sit (4) algebra [92] appears to be 
a good candidate [8| ll 16|.rrT7] for studying the ladder compounds, since its 
Hamiltonian is that of the standard Heisenberg ladder with an additional bi- 
quadratic spin interaction. The physical interpretation of the biquadratic spin 
interaction term is of various kinds. In the continuum limit it represents an 
interchain coupling of spin-dimerization fields which therefore influence the 
low-energy properties of the spin ladder [91]. In the strong coupling limit, the 
effective biquadratic spin-spin coupling arises from virtual electron transitions 
between neighbouring magnetic atoms. This coupling is in general very small 
in spin chain systems [118] , smaller than for example the quadratic interaction 
term. In principle, the strength of the biquadratic interaction in ladder materi- 
als might be accurately determined through calculation of exchange integrals. 

For strong rung coupling, i.e., 4Jy <C J±, the spin-spin exchange interaction 
is dominated by the strong rung interaction. It follows that the contribution 
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from the interaction along the ladder legs, including the biquadratic term, to 
the low-temperature physics is much less than the contribution from the rung 
interaction. However, the leg interactions and the biquadratic spin interactions 
do cause a shift in the critical field values. This kind of shift in the critical 
field values can be balanced by adjusting the interaction strengths. We shall see 
clearly that the values of the coupling constants Jii for the various compounds 
determined from the integrable model are always smaller than the values fixed 
by the standard Heisenberg spin ladder, as a result of including the biquadratic 
interaction term. The integrable spin ladder model provides a consistent and 
elegant way to obtain the ground state and thermodynamic properties of the 
strongly coupled ladder compounds, i.e., with 4Jy <C J±. We emphasise again 
that we set aside any issues regarding the reality of the biquadratic interaction 
term. The advantage of the integrable model approach is that it opens up 
the exact calculation of finite temperature thermodynamic properties via the 
well established methods developed for integrable systems [11[5], including the 
more recently developed high-temperature expansion method reviewed in this 
article. 

The interaction Jn plays an essential role in opening up a Luttinger liquid 
magnetic phase between the singlet and fully-polarized phases for H c \ < H < 
H C 2. Although the physical properties are not sensitive to the value of Jn at 
large temperatures, Jii has a significant influence on the low temperature be- 
haviour. It is argued that for the integrable model the collective effect from 
the biquadratic interaction and the exchange interaction along the ladder legs 
are equivalent to the pure Heisenberg type interaction in the standard Heisen- 
berg spin ladder model for the strong coupling region 4Jn <C J±. The coupling 
constants can be fixed by an overall fit of the critcal fields and thermodynamic 
properties with the experimental results. In this way the integrable model is 
seen to exhibit similar critical behaviour, i.e., in the same universality class, 
in the correlation functions, magnetization, susceptibility, specific heat etc to 
the non-integrable Heisenberg ladder, which has so far been used to examine 
spin ladder materials [20] . 

The Thermodynamic Bethe Ansatz (TBA) |119U120j has been widely applied 
to many physical problems, e.g., to the Heisenberg magnet [12011121] . the 
Hubbard model [EllCm], the super symmetric t-J model [TM[T25l[T26l[T27] 
and the Kondo problem [128] . It has most recently evolved into the exact 
High Temperature Expansion (HTE) method [H^[TB0l[T3Tl[T32l[T33] . The key 
ingredients of the HTE method are the Quantum Transfer Matrix (QTM) 
plli:-i5yi.S6pi.S7Ui:-i8yi.S9U14()pi4iyi42lfTi3] and a functional relation called the 
T-system |144tll45|ll46j . from which one derives nonlinear integral equations 
which can be solved in an exact perturbative fashion. To date this approach 
has been applied to the Heisenberg model [129] . the osp(l\2s) model |130j . the 
su(N) Uimin-Lai-Sutherland model [131] . the higher spin Heisenberg model 
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j!47j . the su(N) Perk-Schultz model [148] and the su(m\n) Perk-Schultz model 
[149] . In particular, in this way, i.e., via the HTE expansion, it has been 
demonstrated that the integrable su(4) ladder model can be used to study 
the thermodynamics and magnetic properties of the strongly coupled ladder 
compounds [8]. 

For strongly coupled ladders, the rung coupling J± can overwhelm the in- 
trachain interaction Ju such that a pair of spin-5 magnetic ions along each 
rung forms an effective spin-2,S ion, allowing the strongly coupled ladders 
to be mapped to one-dimensional spin chains. The integrable ladder models 
based on su(N) symmetry provide a tunable way to describe experimental 
compounds for which the rung coupling dominates the exchange interaction 
along the legs. In a certain basis, the rung interaction terms are diagonal while 
the basis maintains the su(N) permutator exchange interaction along the legs. 
More precisely, the integrable ladder models based on su(N) symmetry can 
be mapped to one-dimensional su(N) [HJIMlEIllSEllHH] chains with additional 
chemical potentials. Consequently, the thermal and magnetic properties can 
be systematically derived from the exact TBA and HTE approaches. 

The TBA method is seen to be most convenient for predicting critical fields 
of relevance to quantum phase transitions in ladder compounds |116yil7|ll50^ 
I151|,ll52j . spin chains [15311154] with strong single-ion anisotropics |155|.ll56j 
and spin-orbital models [15711158] . It has been shown [1161I117|I151] that in 
the strong coupling region the integrable two- leg spin-^ ladder model exhibits 
three quantum phases: 

(i) a gapped phase in the regime H < H c \ = J± — 4 Jn , 

(ii) a fully polarized phase for H > H C 2 = J± + 4Jm, 

(iii) a Luttinger liquid magnetic phase in the regime H c \ < H < H C 2- 

The critical behaviour in the vicinity of the critical points H c \ and H C 2 is 
seen to be of the Pokrovsky-Talapov type. On the other hand, the HTE 
method gives the temperature-dependent free energy, from which physical 
properties such as the magnetic susceptibility and the specific heat can be 
derived [8l ll52j . Indeed, as we shall discuss here, in certain regimes, the 
integrable two-leg ladder with strong rung coupling [8, 152] describes some 
experimental compounds, e.g., (5IAP) 2 CuBr 4 -2H 2 [12], Cu2(C 5 Hi2N 2 )2Cl4 
HUES], (CsH^N^CuB^ [37]) as well as the double spin chain systems 
KCuCl 3 [?,?S[5l] TlCuCl 3 [501E3E81E0], organic polyradical ladders BIP- 
BNO [ISnKIlO], [Cu2(C 2 02)(CioH 8 N2)2)](N0 3 )2 [EE] and the ferrimagnetic 
mixed spin ladder PNNBNO [TIlEHllTJ . 
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1.4 Key results and outline 

In this article we review the calculation of the physical properties of an inte- 
grable two- leg spin-, ladder and an integrable mixed spin- (|,1) ladder. For 
the integrable spin-| ladder with strong rung coupling the singlet ground state 
is separated from the lowest magnon (triplet) excitation by a finite energy gap. 
Notably, the values for the energy gap, A = J± — 4Jji , and the critical points, 
H c i = Jj_ — 4 J|| and H C 2 = J± + 4 Jn , derived from the TBA are in ex- 
cellent agreement with the available experimental results for strong coupling 
ladder compounds. The compounds discussed are (5IAP) 2 CuBr4-2H 2 [32], 
Cu 2 (C5Hi2N 2 )2Cl4 [441145] . (C 5 Hi2N) 2 CuBr4 [47], the organic polyradical lad- 
ders BIP-BNO [T5MT60] and [Cu 2 (C 2 2 )(CioH 8 N2)2)](N0 3 )2 HSJJ. The finite 
temperature thermal and magnetic properties are examined for these com- 
pounds via the exact HTE, with good overall agreement between the the- 
oretical predictions and the experimentally measured thermal and magnetic 
properties. 

For the integrable mixed spin-(i, 1) ladder, the ground state lies in a gap- 
less phase and a one-third saturation magnetization plateau occurs as the rung 
coupling exceeds J± > | Jy In 2. In the absence of a magnetic field the model ex- 
hibits three quantum phases associated with doublets, quadruplets and su(6) 
symmetry. The fractional plateau corresponding to a fully-polarized doublet 
state opens at the critical field H c \ and vanishes at the critical field H C 2. If 
the magnetic field is greater than H c %, the quadruplets are fully polarized. 
The phase diagram is reminiscent of the spin-| chain |153|]154] . The magnetic 
properties of the organic mixed ferrimagnetic ladder compound PNNBNO are 
examined via the HTE. 

The ground state properties of two integrable spin-orbital models are also 
treated via the TBA method, from which it is seen that the single- ion 
anisotropy and orbit splitting field can trigger a magnetization plateau with 
respect to different Lande factors with spin and orbital degrees of freedom 
[157^158] . The magnetization for spin and orbital degrees of freedom is also 
calculated numerically from the TBA equations. 

The article is organised as follows. In section [2J we review the two inte- 
grable spin ladder models - a two-leg spin-^ ladder and a mixed spin-(i, 1) 
ladder - along with their Bethe Ansatz solutions. In section [3l we discuss 
the TBA equations for the two- leg spin-^ ladder and the investigation of the 
ground state properties via the TBA equations at zero temperature. The ex- 
act and resulting numerical calculations for the magnetization and suscepti- 
bility are treated in detail. The Quantum Transfer Matrix method is reviewed 
in section HJ along with the derivation of the nonlinear integral equations 
in terms of the T-system. It is then shown how the free energy is derived 
from the HTE. The results are then applied in section [5] to examine the 
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thermal and magnetic properties of the strong coupling ladder compounds 
(5IAP) 2 CuBr 4 -2H 2 0, Cu 2 (C 5 Hi2N2) 2 Cl4, (C 5 H 12 N)2CuBr4, and the organic 
polyradical ladders BIP-BNO and [Cu 2 (C 2 02)(CioH 8 N 2 ) 2 )](N03) 2 . In section 
[fH the same approach is applied to the mixed spin-(^, 1) ladder with a discus- 
sion of the ladder compound PNNBNO. Section [7] is devoted to the study of 
the magnetic properties of the spin-orbital model. A concluding discussion is 
given in section [8j 

It should be noted that readers more interested in the physical properties 
of the ladder compounds, rather than their detailed calculation, may choose 
to skip sections [3J and [H In any case the relevant results are collected at 
the beginning of section [5] on the comparison of the thermal and magnetic 
properties with the experimental data. Conversely, readers more interested in 
the mathematical methods may prefer to concentrate on sections [3J and [H 



2 Integrable spin ladder models 



The general two-leg ladder configuration is depicted in figure [2j The two-leg 
Heisenberg spin-i ladder model is defined by the Hamiltonian [T8 1 I20 1 I63] 



L L L 

n = J \\ E$ • ^i+i + f i • + J± E ^ • ^ " ^(Sl + Tf). (2) 
j=i j=i j=i 



Here Sj = and Tj = (Tj,Tj,T?) are spin-i operators acting on 

site j. For example, Sj = where aj = (cr?, crj, cj|) are Pauli matrices. The 
essential parameters are the rung coupling J±, the intrachain coupling Jy and 
the external magnetic field H. The Bohr magneton is fi-Q and g = g s = gt is the 
Lande factor. There are L rungs and periodic boundary conditions, Sl+i = Si 
and Tl+i = Ti, are imposed. 




Figure 2. The two-leg spin ladder. In general, sites on the lower (upper) leg carry spin-5 (T) with 
g s (gt) the Lande factor along the leg. Jy is the intrachain or leg coupling and J± is the interchain 

or rung coupling. 



This model has been extensively studied via numerical methods [T8 jll62l 



12 



M.T. Batchelor, X.-W. Guan, N. Oelkers and Z. Tsuboi 



I163yi644 l27]. perturbative field theory [S3"1IM1 I(37 ]I165P68| and transfer matrix 
methods [66]. Some essential physics of the model can be obtained from the 
strong coupling region J± S> Jii . Consider the more general spin-S 1 ladder. 
Then in this limit the ladder effectively becomes a spin-25 chain with single- 
ion anisotropy. The physics of the ladder /chain equivalence holds even for 
weak coupling, but this is by no means obvious. The gapped phase persists for 
arbitrary rung coupling, with the spin-spin correlation decaying exponentially. 
The analogy between the physics of the spin chain with 2S even/odd and the 
spin ladder with even/odd number of legs is well described by the quantum 
field theoretic approach |166yi67|ll68| . Both the antiferromagnetic Heisenberg 
spin chain and the ladder can be mapped to the nonlinear sigma models. 
The field theory approach thus gives a unified description of the Heisenberg 
chains and ladders. In this framework the gapped and gapless phases of the 
spin chains and ladders are governed by the topological term 9 = or it, 
respectively, in the nonlinear sigma model [169, 168J . 

2.1 The integrable two-leg spin-^ ladder 

The simplest integrable spin ladder is defined by the Hamiltonian [92] 

L L 

H = J,, Wfcg + Jl % ■ Tj - VbqH YJ,SI + Tf) (3) 

3=1 3=1 

where 

L 

Wicg = E • Si+i + Tj ■ f j+1 + 4(5, • S j+1 )(fj ■ f i+ i)) . (4) 

3=1 

Periodic boundary conditions, Sl+i = Si, Tl+i = Ji, are applied. This inte- 
grable Hamiltonian contains a biquadratic term, which minimises the Haldane 
phase [9T |lll6pil7j and causes a shift of the critical value of the rung coupling 
J_i_ at which the gap vanishes (recall that the gap vanishes at zero rung cou- 
pling for the Heisenberg ladder). The other key feature of the model is that 
it has the same magnetic field and rung interaction terms as the standard 
Heisenberg ladder ([2]). In the absence of symmetry breaking potentials, the 
dispersion of the low-lying excitations is entirely caused by the leg interaction. 
However, symmetry breaking potentials, such as rung interaction or external 
field, result in different band fillings such that the lowest spin excitation is 
gapped. The two-body and biquadratic interactions along the legs are effec- 
tively quenched by the strong rung singlet. The model thus lies in the same 
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phase as the standard Heisenberg ladder. 

The critical behaviour of the model is determined by the competition be- 
tween the rung and leg couplings and the magnetic field H. A key factor in 
the integrability of this Hamiltonian is that 7i\ cg is (up to a constant) simply 
the permutation operator corresponding to the su(4) algebra [11441115] . This 
can be seen by performing a change of basis to the states 



|i> = ^(IU>-|it»,|2) = |TT> 
|3> = ^(IU) + UT)),|4) = |u> 



(5) 



where state |1) denotes the rung singlet and the remaining three states denote 
the triplet components. Here we have set | f) = (J) and | j) = Q). Thus 

ITXTI = 1(1 + **), = 1(1-0, ITXII = <? + and IIXTI = a'. Note 

that in this new basis the leg part remains of the same form while the rung 
term becomes diagonal. The rung term reduces the su(4) symmetry of TC\ cg 
to su{3) © it(l) symmetry under the action of the local chemical potential 
J±\l)jjiM- Switching on the magnetic field further breaks this symmetry due 
to Zeeman splitting. Moreover, the su(3) © u(l) symmetry does not survive if 
we consider different g factors along the two legs. Let gt and g s be these Lande 
factors along the upper and lower legs, respectively. Then the Hamiltonian 
reads fT571fT58l 



L L 

H = J||Hi cg + Jl ]T Sj ■ fj - vbH Y,(9sS* + g t T?) (6) 
i=i j=i 

where the local eigenbasis on the jth rung is given by 

^ = 7r =^(lU)^ ( - ) liT)) ) <h- = I TT) 

™(lU)+ y (+) IIT)), 3 = |U) (7) 



with 



= (g s - gt)» B H/J ± ±^Jl + (g s - g t f ^ j J\. (8) 

In this procedure it is important to note that the leg content, Hi eg , equation 
(j3|), of the Hamiltonian is invariant under any change of rung eigenbasis. It 
is also invariant under any choice of reference state for the application of the 
algebraic Bethe ansatz, while the rung and the field terms are altered by these 
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changes. The validity for reordering the basis arises from the fact that the 
intrachain interaction is the permutation operator, which is independent of 
the choice of basis for the vector space V = V\ © V2. On the other hand, the 
rung interaction and magnetic field vary the energy levels of the eigenbasis 
states following from equations © and ([7]). 

Generally speaking, the difference between gt and g s is often very small and 
even negligible for the isotropic two-leg ladder model, but this may be different 
for the spin-obital models |170pi7ipi72j . For simplicity, unless otherwise stated, 
we consider here the isotropic ladder model ^ for the application of the TBA 
and HTE methods. 

In section [Hit is shown that the integrable Hamiltonian ([3]) can be derived, 
up to an unimportant constant, by exploiting the spectrum parameter v and 
the row-to-row transfer matrix t(v) = tro-Ro,L(^)-Ro,L-i(^) ■ • • Ro,2(v)Ro,i(v) 
via the relation 

, L L 

H = J,, — ln*(«)| w= o- J±^|l) i (l|,--// B ^E(|2) i (2|,-|4) i (4| i ) (9) 

in terms of the underlying su(4) i?-matrix, which is simply R(v) = vl + V, 
where V is the 16 x 16 permutation operator (cf. equation (|79p ). After some 
algebra, the energy spectrum can be obtained from relation Q as |116|I117] 

Mi 1 

^ = - J llE^n - J±N Q - » B g(N+ - N_)H (10) 

i=i (vl >) 2 + 3 

in which iVo denotes the number of rung singlets and N + and iV_ denotes the 
number of the triplet components |2) and |4), respectively. 
The number M\ appearing in the energy spectrum (|10p is the number of 

Bethe roots vj- which satisfy the r = 3 case of the general su(r + 1) Bethe 
equations 

M k -! (k) (fc-1) ,i- M k (fc) (fc) , ■ M k +i (fc) (fc+1) 1- 
V) — V: + ^1 t-4- V) — Vf +1 V) — V, — -SI 



n j. * ^_ = tt-^ - — tt - 1 - - (in 

11 Jk) Jk-1) 1. 11 Jk) n ,( fc ) ; Al n ,( fc ) „,( fc+1 ) , 1; 

1=1 "' 



i=i v j ~ v i -2 1 1=1 v j ~ v i "i 1=1 v j ~ v i +1 [ 



obtained by application of the nested Bethe Ansatz [2], [9j [TOl HT] . Here k 
1,2, ... ,r and j 
and M = L. 



1,2, ... , r and j = 1, . . . , M^, with the conventions = v^ r+1 ^ = 0, M r+ i 
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If g s / gt, the corresponding energy is given by [157^158] 



(12) 

with h = /j,bH/J±. In this case N^ a (a 6 {0,1,2,3}) is the number of (f) a 
states. 

2.2 The integrable mixed spin-(^,l) ladder 

Recalling figureEl the integrable mixed spin-(i, 1) ladder model is defined by 
the Hamiltonian [TTilfTTHlfTHn] 

L 

-VBHY.igtTf+gsS*) (13) 

1 + Sj ■ S j+1 + (Sj ■ S^) 2 ) . (14) 

For this model Tj remain spin-^ operators but now Sj are spin-1 operators 
and we consider general Lande factors gt and g s . Periodic boundary conditions 
are again imposed and L is the number of rungs. 

By construction 7i\ e g defined in (|14p is the permutation operator with su(6) 
algebraic symmetry. In this case the rung term in (113f) breaks the su(6) sym- 
metry into an su(4)(Bsu(2) symmetry, i.e., quadruplet plus doublet symmetry. 
Now changing the fundamental basis of su(6) into rung quadruplet and dou- 
blet states (via Clebsch-Gordan decomposition) results in the six-dimensional 
state space splitting into the direct sum of quadruplets and doublets with re- 
gard to the rung interaction. The projectors onto the doublet and quadruplet 
subspace are given by 

V d = -l(f -5-i), 7> q = §(f -S + l). (15) 

Although the magnetic field preserves the integrability of the leg part of the 
Hamiltonian, the different ^-factors on each leg break the doublet/quadruplet 
basis for the Hamiltonian (fl~3j) due to Zeeman splitting. Fortunately, there is 



S 



M 1 

-'»£■ 



i>f } ) 2 + i 



5 J± [l + y/l + (9s 



9t 



\ 2 h 2 



\M9s + gt)HN^ + ±ti B (9s + 9t)HN^ 
\J L \\ - Vl + (9s ~ 9t) 2 hA 



L 

Tt = J||?Yicg + J± ^ ^ Tj • Sj — 

L 
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another basis, 



, 1 (|i,-i) + i/P°|Q,j)) 

. /1 j. f„(±h2 V 2 7 



i + (vl 

2 

^ = I 1 (l - 1, |) + 2/SlO, -I)) (16) 



t/>3 = |1, \ 



i 2i 

V>-2 = |-1,-|> (17) 



that diagonalises the rung and magnetization terms simultaneously. Here the 
quantities are given by 



yj& = av^[( 5 s - 9t)h' + a] ± ^ 1 + - + a) 2 (18) 

where a = dr, and // = /^ibH/J±. Further 




11) = I I , |0) = I 1 I , I - 1) = I I , |i) = (l) , I - \) = (J) (19) 



with \b,a) = \b) (g) |a). If g s = gt, the basis states tp[ and V'ii reduce to the 

2 2 

doublet, with the other states reducing to the quadruplet. It follows that the 
leg and rung parts of the Hamiltonian (I14j) are given by the relation 

H = J\\ — ]nt(v)\ v=0 +^^^ Q) |^ Q) )i(^ Q) Ij + constant (20) 

i=i 

associated with the row-to-row transfer matrix. Here a = +,0,— and (3 = 
±|, with tfig 1 ^ = ipp. The chemical potentials /A a ' are given in (|22p below. 
The eigenspectrum is obtained from relation (I20p and the Bethe Ansatz to 

be nsg 

f = - J iiE^^-+^ ( 21 ) 
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where the Bethe roots satisfy the Bethe equations (jlip with r = 5. The 
energy contribution from the rung interaction and the magnetic field terms is 
given by 



£±+h 



\J± ~ \9s^bH - ^J±\/l + \{g s h' - g t ti + §) 



1\2 



(-) 



+ 



\J± + \9sVbH - -^J±\f l + 5 (ft.fr' - g t ti 



1\2 



+ [\JL-{\gt+9s)^BH] N, 



+ 



+ 



ijjL - ^.MB-ff + 75 J ±V 1 + Udsti - g t ti + i) 2 



(+) 



1^1 + hs^sH + -^J±\/ 1 + 5 - fftfr' 



1\2 



+ ^J± + (^t + ft)AiB^] 



(22) 



In the above equation, iv£^ and are the occupation numbers of the 

2 2 

corresponding states and /?/ = /ibH/J±. 



3 Thermodynamic Bethe Ansatz 
3.1 Derivation of the TBA equations 

In this section we turn to the study of the ground state properties and the 
phase diagram of the integrable two- leg spin-^ ladder ([3]) by means of the 
Thermodynamic Bethe Ansatz (TBA). In the thermodynamic (L — > oo) limit, 
the roots of the Bethe equations for the ground state are real. At zero tem- 
perature there are no holes in the ground state root distribution. Yang and 
Yang |119j introduced the TBA in treating the thermodynamics of the one- 
dimensional boson model with delta function interaction, where the elemen- 
tary excitations involve moving quasimomenta out of the Dirac sea. In general 
complex conjugate pairs of Bethe roots form bound states of two magnons. 
The string hypothesis has been developed to classify the many permissible 
bound state solutions to the Bethe equations, originally for the Heisenberg 
chain |120j . the one-dimensional Hubbard model [122] , the supersymmetric 
t - J model P21H251H2E1II2Z!. Kondo problems p8] and the higher-spin 
Heisenberg chain |173] among others. Recall that the Bethe equations (llip for 
the ladder model (J3]) consist of a set of three coupled equations with Bethe 

(k) 

Roots of three different colours, , k = 1,2,3. Taking the thermodynamic 
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limit |119yi2U|ll74j results in the string solutions 



v (k) n . = y (k) n 



_ + |i(n + 1 - 2j) (23) 



a k j " a k ' 2 

where j = 1, . . . , n and a& = 1, . . . , N n k \ Here u^'" , A; = 1, 2, 3, are the 
positions on the real axis of the centre of an n-string, where n is the string 

(k) 

length in each colour. There are N n strings of the n-string form. The number 
of n-string satisfy the relation = YlnLi 71 -^^, where is the 

total number of Bethe roots of colour k. 

It is assumed that the distributions of Bethe roots are dense along the real 
axis. The missing roots of each colour k are the holes, with positions v^ h a 

on the real axis. The distribution densities of n-string roots, pn (v), and holes, 
p^ h (v), are defined to be 

P W(V) = Km (k)n 1 nXnT (24) 

pW h (v) = lim - — ^ 1 — tj- (25) 

L^L( v m n ak+ i-v {k K k ) 

for the three colours k = 1,2,3. 
Adopting the standard notation |173| 

A nm (\) = 5{X)5 nm + (1 (A) + (A) 

Min(n,m) — 1 

+ 2 ^ a \n-m\+2lW (26) 
(=1 

Min(n,m) 

Onm(A) = E a n+rn+ i_2«(A) (27) 
Z=l 

with a n (A) = 2^ n 2y4 + ^2 , the Bethe equations (fTTjl can be transformed after a 
long though standard calculation into the form 

p^\v) = a n (v) - J2(A nm * p$)(v) + J2(*nm * (®)(v) (28) 

m m 

/4 2)/ » = - E(^- * + * + (28) 

m m 
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where * denotes the convolution integral (/ * g)(X) = ff° /(A — X')g(X')dX' . 

In order to find the equilibrium state of the system at a fixed temperature T 
and external magnetic field, the free energy F = E — TS — p^,gM z is minimised 
with respect to the densities. To this end, the energy per site can be written 

as 

°° /'OO 

S/L = - J,, ]T / dva n (v)p£\v) - J ± N (31) 

n=l^-°° 

in terms of the densities of roots and holes. Following |119j the total entropy 
of the system is given by 

00 poo I 3 

S = E / dV \ J2(Pn ] («) + Pn k)H (v)) Hp ( n k) (v) + P^ "(«)) 
n=1 J-oo \ k=l 

-$\v)]n$\v)-pW\v)]npP h (v)) (32) 
with magnetization 

M z = fj, B g(N + - jV_). (33) 

Minimizing the free energy, i.e., 5F = 0, and making use of the equations 
(|2"S]) -([50 ]l gives the TBA equations 

^ 00 
ln(l + v W(v)) = -gi h \v) + A nm * ln(l + 

m=l 

00 

-^a nm *ln(l + i fc - 1 i" 1 )(i;) 

m=l 

00 

-^a nm *ln(l + ^ +1 )- 1 )( W ) (34) 

m=l 

in terms of the functions rjn\v) = p^ h {v) / p n k \v) . Henceforth for nota- 

tional simplicity we will omit the argument v. The driving terms gn\v) = 
nx ( k > — 5i^27ra n (v) appearing in the TBA equations ([34"]) play a key role to 
be discussed further below. Here is the chemical potential for the colour 
/c-branch. 

The TBA equations are often written in terms of the dressed energies e n k \v), 
defined via r£?\v) = exp(ei fc) (v)/T) for k = 1,2,3 with {v) = e { n ] (v) = 0. 
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OO 



Specifically, 

e[ k) = g[ k) + Ta 2 * ln(l + e"^^) + T(a + a 2 ) a m * ln(l + e ~ e ™V T ) 

m=l 

oo 

-T *m * (ln(l + e-^~ l) / T ) + ln(l + ' T )) (35) 



m=l 



with 

e (f) = 5 « + Tax * ln(l + e- e »-/ T + Ta 2 * ln(l + e^'^ 

oo 

+T(a + a 2 ) a m _ n * ln(l + eT e ™ l T ) 

m=n 

oo 

-T £ a m _ n+1 * (ln(l + e-^ 15 ^) + ln(l + e"^/ 3 ")) (36) 

m=n 

for n > 1. The dressed energies play the role of excitation energies measured 
from the Fermi energy level for each colour, and satisfy a Fermi-Dirac-like dis- 
tribution. This can be seen from the ratio of the number of occupied vacancies 
to the total number of vacancies in the interval v to v + dv for colour k. Namely 

A%) _ i (37) 



i +«#'<<•>' 

The dressed energies tn\v) for k = 1,2,3 are monotonically increasing func- 
tions of \v\ and even, i.e., e$(v) = en (—v) with 

Urn f=-W=2x<*>. 

n— >oo fl 

In the above equations the driving terms gr^ and depend on the choice 
of the basis order due to the energy level splitting among the singlet and the 
triplet components with respect to the competition between the rung interac- 
tion and the magnetic field, see (|10p and (|12p . It is worth noting that for the 
ladder Hamiltonian ([3]) the singlet rung state is energetically favoured when 
J_l_ > (antiferromagnetic coupling), while for J± < (ferromagnetic cou- 
pling) the triplet rung state is energetically more favourable. Here we consider 
only the antiferromagnetic coupling regime, i.e., with Ju > and J± > 0. Al- 
though the physics of the model doesn't depend on the basis order, the driving 
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terms vary the Fermi surfaces of each branch in different ways. It is most con- 
venient to order the basis according to the energy levels of the singlet and 
triplet components which provides a natural order with respect to the band 
filling. For strong rung coupling there are only two components, the singlet 
and the lowest component of the triplet, which compete in the ground state. 
The other components induce a gap above the ground state. In this natural 
basis order, the TBA equations for the ground state reduce to one band with 
finite Fermi points. 

Explicitly, in the strong coupling regime, i.e., J± > 4Jn and J± > H, the 
basis order is chosen as (|1), |2), |3), |4)), corresponding to the Bethe states with 
number N\ of state |1), up to number N4 of state |4). Here N\ = L — M±, N2 = 
Mi - M 2 , N 3 = M 2 - M 3 and N4 = M 3 (see equation (pi])). The driving terms 
are thus given by 

g{ 1] = -J||27roi + J± ~ HBgH, gf ] = gf = (i B gH. (38) 
The higher order driving terms are gn^ = n^-QgH, gffl = ng^ and g n ^ = 

(3) 

ng\ . Subsequently, up to a constant, the free energy is 

/oo 00 
;>> n (t;)ln(l + e- e » ^' T )dv. (39) 
-°°n=l 

However, if H > J±, the basis order can be taken to be (|2), |1), |3), |4)) for a 
simpler analysis of the TBA equations. The driving terms are then 

g^ = -J||27rai- J ± +f, B gH, g® = J ± , gf ] = mg H (40) 

with the higher order driving terms remaining unchanged. In this case the free 
energy is 

/oo 00 
]Ta» ln(l + e- £ " W T )dv. (41) 
-°°n=l 



3.2 Ground state properties and phase diagram 

We turn now to the ground state properties and the phase diagram of the 
integrable ladder model, as revealed by the analysis of the TBA equations. 
In the low-temperature limit T — > only the negative part of the dressed 
energies e®, denoted by = min{e^,0}, contribute to the ground-state 
energy. Here we use the simplified notation = ef . The TBA equations 
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then become 







— 02 * 


+ ai * e (2) 




e (2) 




- a 2 * e^ 2) ~ 


+ ai * 


e (l)- +e (3)-" 


(42) 


£ <3) 




- a 2 * 


+ ai * e (2) 





where <?[ , a = 1, 2, 3, are the driving terms with respect to the different basis 
orders given in (f38j) and ((101) . 

In the absence of the symmetry-breaking chemical potential and magnetic 
field terms, the Fermi seas of the three bands are completely filled and have 
no Fermi surface, as depicted in upper left panel of figure [3l The Fermi sea of 
each band may be raised or lowered by the rung interaction and the external 
magnetic field, with spinon excitations created by flipping a spin. The chemical 
potential, rung interaction and magnetic field terms vary the Fermi surfaces 
at the bands, leading to interference effects between the two Fermi points. 

In the antiferromagnetic regime ( J± > 0) it follows that if H = the triplet 
excitation is massive, with the gap given by A = Jj_ — 4 Jn . This configuration 
is depicted graphically in the upper right panel of figure [31 where all of the 
dressed energies are gapful. The critical point J+ = 4Jn indicates a quantum 
phase transition from the three branches of the Luttinger liquid phase to the 
dimerised u(l) phase. The magnetic field lifts the fermi seas of and e^ 3 ), 
which can be seen from the TBA equations (1421) . 

For Jj_ > , it can be shown that the two components of the triplet states, 
|3) and |4), do not enter into the groundstate for high magnetic field. However, 
the magnetic field shifts the energy level of the triplet component of |2) closer 
to the singlet groundstate. The gap can be deduced via the magnetic field H. 
The first critical field occurs at H c i, where g[iBH c \ = A, i.e., the effect of 
the magnetic field is to close the gap. The quantum phase transition is from 
a gapped to a gapless Luttinger phase. When H > H c \ the triplet component 
of 1 2) begins to enter into the ground state with a finite susceptibility. The 
dressed energy configuration for this case is depicted in the lower left panel 
of figure (3j The dressed energies e^ 2 ) and e^ 3 ) are gapful, whereas e^ 1 ) has 
two Fermi points. The excitation energies are proportional to the momentum 
of particles. If H > Hip = J_l//Ub<7, the triplet component of |2) becomes 
the lowest energy level with driving terms (j40p . At the point H = -ffip, the 
magnetic field completely diminishes the rung interaction so that the Fermi 
sea of the dressed energy e^ 1 ) is completely filled, whereas the others are still 
gapful (see lower right panel of figure [3]) . In this way we see that the ground 
state is in a fully-polarized ferromagnetic state when the magnetic field is 
greater than H C 2 = J± + 4Jm . 
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Figure 3. Depiction of the dressed energies for the su(4) TBA equations [142 1 . The four panels 
indicate the quantum phase transitions which occur in the ground state of the integrable ladder 
model (J3j with strong rung coupling. The horizontal axis is the spectral parameter with the 
origin ( v (*0 = 0) at the minima. The upper left panel shows the dressed energies in the absence of 
the symmetry-breaking chemical potential. In the presence of the chemical potential, the basis is 
reordered so that the Fermi surfaces of the dressed energies remain in the order e^ 1 ) lowest, 
second lowest and e' 3 ' highest (see text). 

The TBA predictions for the occurence of critical fields H c \ and H C 2, along 
with their values, are in good agreement with the experimental results for sev- 
eral strong coupling ladder compounds [11611117] . as will be discussed further 
in sections 13.31 and [5j 

3.3 Analytical and numerical solutions of the TBA equations for the 
magnetization 

The TBA equations gSJ) may be solved numerically [11601170157015801520150] 
to study the magnetization of the integrable spin ladder model in detail. As a 
concrete example we take the strong coupling compound (5IAP)2CuBr4 -2^0 
[42] for which the coupling constants can be used in the numerical calculations. 
This compound was recognized as a strong coupling ladder compound with the 
critical field values H c \ ~ 8.3T and H C 2 ~ 10. 4T. For the integrable spin lad- 
der model ([3]) the critical field values are H c \ = J± — 4 Jn and = J± + 4 Jy . 
In this way reasonable values for the coupling constants Jy and J± can be sug- 
gested from the experimental critical field values. The g-factors are uniquely 
determined by experiment and are not free parameters in this approach. How- 
ever, the coupling constants fixed by the critical fields only lead to a good fit 
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with the energy gap. In order to make a good overall fit for each of the key 
thermal and magnetic properties, the coupling constants should be fixed by 
the energy gap, the magnetization and the susceptibility. Once the coupling 
constants are fixed, all thermal and magnetic properties follow naturally. For 
this compound a fit to the experimental critical fields and the magnetization 
at different temperatures suggests the coupling constants J± = 13.3 K and 
Jii = 0.2875 K, where we have used g = 2.1. The leg coupling constant J» is 
smaller than for the pure Heisenberg model leg interaction due to the pres- 
ence of the biquadratic term, which strengthens the overall leg interaction. 
Recall that the leg interaction ^ differs from that for the pure Heisenberg 
ladder model. As a result the leg coupling constant Jy for compounds deter- 
mined from the integrable model (j3]) differs from that determined from the 
Heisenberg ladder model ([2]). In order to compare the difference between Jy 
determined by these two theoretical models, a scaling parameter 7 was intro- 
duced [11611117] . where J\\/"f characterizes the interaction strength along the 
legs for the integrable spin ladder. 

From the above analysis, when H < H& = 8.6 T the ground state is the 
product of rung singlets such that the magnetization M z = 0. If H& < H < 
Hip = 9.42 T, the dressed energies and e^ 3 ^ are both gapful, so that the 
magnetization is given by 



M z = / p^(v)dv (43) 



where 



2irv 2 + \ 2ir y_ Q 1 + (v - k) z 
In the above equations Q denotes the Fermi boundary determined from 

eW(u) = -2irJ«ai(v) + J± - fi B gH - / — e^~(k)dk. (45) 
11 2vr J_q l + (v- k) 2 

In the region where the magnetic field is greater than 9.42 T but less than 
H C 2 = 10.2 T, the Fermi boundary is determined by 

e«(?;) = -2nJ liai (v) - J ± + fi B gH - ±- [ Q — , 2 (k)dk (46) 

2vr J_q l + {v - ky 
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Figure 4. Comparison between theory and experiment for magnetization versus magnetic field H 

in units of saturation magnetization for the strong coupling ladder compound 
(5IAP)2CuBr4 • 2H2O |42j . The theoretical curves are obtained from the integrable model with the 
values Jj_ = 13.3 K, J11 = 0.2875 K and g = 2.1. The full curve is obtained from the TBA equations 
at zero temperature. For further comparison, the dashed curve is obtained using the HTE at 
T = 0.4 K, as outlined in section [5] The experimental magnetization at T = 0.4 K (diamonds) is 

shown for comparison. The discrepancy in the curves at T = 0.4 K is due to paramagnetic 
impurities, which result in nonzero magnetization in the singlet ground state. At the inflection 
point, H = Jj_/Mb<?, t ne magnetization is 0.5. Both the TBA and HTE curves indicate values of 
the critical fields in excellent agreement with the experimental results. 



and the magnetization is 

rQ 

M z = 1 — / p {1 \v)dv. (47) 
J-Q 

The numerical solution of these equations gives a reasonable magnetization 
curve (see figure U]) which passes through an inflection point midway between 
H c i and H C 2. The inflection point has a physical meaning in that the proba- 
bilities of the singlet |1) and the triplet state |2) in the ground state are equal. 
At zero temperature, in the strong coupling regime, the one point correlation 
function, {Sj -Tj) = — |, lies in a gapped singlet groundstate, which indicates 
an ordered dimer phase, while (Sj-Tj) = | in the fully-polarized ferromagnetic 
phase. However, in the Luttinger liquid phase, {Sj ■ Tj) = — | + M z . This in- 
dicates that the magnetic field increases the value of the one point correlation 
function. 
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3.4 Small Fermi boundary expansion method in the vicinity of the 
critical fields 

It was shown [11611117] that the magnetization in the vicinity of the criti- 
cal fields H c \ and H C 2 depends on the square root of the field, indicating a 
Pokrovsky-Talapov transition. An analytic expression for the field-dependent 
magnetization in the vicinity of the critical points H c \ and H C 2 can be ob- 
tained by a small Q approximation. The lower left panel of figure [3] shows the 
dressed energies e^ k \ k = 1, 2, 3 in the vicinity of the critical point H c \. 
Performing a small Q approximation in (|45]) gives 

e^(v) = -27rJ,|ai(tO + Jj_ - y, B gH - ±- f 2 e^~(k)dk 

11 2vr J_ q 1 + (v - ky 

--TTT^ + J ^--T^ e(1) (°) ( 48 ) 

1 + 4V Z 7T 1 + V z 

with J e ff = Jj_ — (J-bqH. It follows that 

C (1W ~ ( Jcff I 4J H^ - ^9(H - H cl )n 

l j 7T + 2Q " 7T + 2Q 1 yj 

with H c \ = (J± — AJ\\)/(hb9)- The condition e^\Q) = then gives 



n~ / [Jvg(H - H cl ) 

Q ~\jlQJ\\-h^g{H-H ci y {M) 
Similarly, equation (|44p becomes 

2-k V Z + 1/4 7T(1 + IT J 

with pW(0) = ^rpm ■ The magnetization for small Q can then be obtained as 

W_/V W *«i2(l-f). (52) 

It follows that the magnetization in the vicinity of iT c i has a square root 
field-dependent critical behaviour of the form 

M z n~^B9(H-H cl )/J {] . (53) 
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The thermodynamics may be calculated from the free energy (|39|) and (|41|) . 
To demonstrate this, we look at the free energy in the vicinity of H C 2 via 
the TBA equations (|42p . When H > J±/hb9, the lowest energy level |2) 
dominates other states as the magnetic field increases. The basis order is chosen 
to be (|2), |1), |3), |4)). In fact, a Luttinger liquid phase lies in the regime 
H c i < H < H C 2, where the ground state can be described by an integrable 
XX Z model with an effective field H — J±/hb9, see (|46p . Similarly, in the 
vicinity of H C 2, from (|46p . we have 



/ ^g(H c2 - H) 
Q Y 16J,| - 5^ B g(H c2 - H) 

with H C 2 = (Jj_ + 4J\\)/iib9- Here the Fermi boundary Q is very small for H 
near H. c i- Subsequently, the free energy can be derived from (HTj) as 



/°° (i) 
ai(u)ln(l +e" £l ^' T )dv 
-oo 



oo 



HbqH + / ai{v)e{' (v)dv 



-H B gH+ / ai(v)e\' (v)dv 
J-Q 

rQ 

-fj, B gH+ / ai(v) 
J-Q 



Q 

2irJ\\ ai (v) + J cS - - 1 ; 2 

7T 1 + IT 



= _ wgfl _^(^-H) / 1 + jg\ 

7T V 7T / 

For very small Q the magnetization is given by 
with M s = hb9- Thus 



The results ([33]) and ([35]) demonstrate the square root critical behaviour. 
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They coincide with the numerical results, as applied to the strong coupling 
ladder compound (5IAP)2CuBr4 • 2H2O in figured! 



3.5 Wiener-Hopf method 

The Wiener-Hopf technique has been widely used to solve functional equations 
which are obtained by Fourier transformation of certain integral equations 
|128U124pi25pi75j . In this subsection we demonstrate how to analytically solve 
the TBA equations with a large Fermi surface. 

The higher energy levels of the triplet states 1 3) and 1 4) are always supressed 
in the ground state for strong rung coupling (Jj_ > 4Jy). The spin waves are 
hard-core bosons with massless excitations. In the regime H c \ < H < H C 2, 
the ground state is that of the integrable XX Z chain with an effective field 
H — J±/[j.b9- This field triggers two pure states, namely |0) for H < H c \ 
and |2) for H > H C 2, which effectively correspond to the states ||) and ||) 
of the XX Z model, respectively. The inflection point H = J±/^b9 may be 
considered as the su{2) point because at this point J± is diminished by external 
fields. It reduces to the TBA equation withought any chemical potential for 
XXX model where the ground state is the su(2) Luttinger liquid. These two 
components are degenerate and fully fill the energy band e*- 1 ) . All other bands 
are empty (as per the lower right panel of figure [3]). When the magnetic field 
is greater than Jj_, i.e., when H > the Zeeman energy introduces a 

large Fermi surface of the dressed energy with regard to the basis order 
(| 2), |1), 1 3), 1 4)). Logarithmic corrections are associated with the interference 
effects between the two Fermi boundaries. 

The integral equations (|45[) and (|46p exhibit a large Fermi boundary for H 
in the vicinity of J±/^b9 and can be solved by the Wiener-Hopf method. We 
thus consider the case of H in the vicinity of J±/hb9 such that the ground 
state is described by the su(2) TBA equation with an effective field. With 
J c g = — Jj_ + fiBgH, equation (j36]) can then be rewritten as 

eW( v ) = -J||27rai(» + J cS - («2 * e (1) ")(^)- (59) 

Taking Fourier transforms, defined by f(co) = /(A)e ltjA dA, gives 



eW(o;) = -J||2tt 6 , + \J cS 5{oo) + — ^e (1)+ M- (60) 
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Define the function G(X) by 

G(X) = — I c — - e-' lX "du. 
v ; 2vr 7.^ 2 cosh f 

Then the inverse Fourier transform indicates that the dressed energy satisfies 
a Wiener-Hopf type equation of the form 

-B 

,(1) 



C W(A) 



Because J e g is very small, the Fermi boundary is very large. Let A — > A + -B 
and y(A) = e^(A + 5), then equation (fB"T1) becomes 

r 00 

»< A > - - J ' cosh,(A + B) + W + / °< A + 2B " k ^" k 

pOO 

+ / G{\-k)y{k)dk. (62) 



Now G(A + 2£? — k) — > e _B as -B 3> 1. The driving term can be expanded in 
powers of e _s . We can thus take an expansion y = yi+y 2 +y3 + - ■ ■ with respect 

to the order of e~ mB , m = 1,2, The leading term y\ describes interference 

between the two Fermi points. The energy of the excitation is proportional to 
the momentum. Subsequently, we have 

poo 

|ft(A) = -J|,2vre-( A+ ^ + ±J eff + / G(X - k) Vl {k)dk (63) 

Jo 

poo poo 

y 2 (X) = G(X - k)y 2 {k)dk +/ G(X + 2B - k) yi (k)dk (64) 
Jo Jo 

poo poo 

y 3 (X)= G(X- k)y 3 (k)dk+ G(X + 2B - k)y 2 {k)dk (65) 
Jo Jo 

and so on. 

Each of these equations is of Wiener-Hopf type and can be solved in an ana- 
lytic fashion. Here we content ourselves with the leading term, i.e., considering 
the leading pole in the decomposed functions. The remaining higher terms y n 
can be obtained by iteration. Introduce the functions y±(A) by 

^ If' a So, = (") 
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Taking the Fourier decomposition, we have yi(u) = yf{uj) + y^(u), where 
y 1 (w) are analytic in the upper and lower halves of the complex plane, re- 
spectively. In order to solve the Wiener-Hopf equations, we need to find a 
decomposition of the kernel and the driving term in (163 j) . To this end, we 
denote the driving term by 

<t>(\) = -J||2vre-( A+i3 ^ + ij cff (67) 

and define Fourier transform for driving term 



<MA) 



0, A < 

£ f e-^+^du = </>(\), A>0 

rh m = J ' A>0 

\&fe- iuX 4-(u;)<kj = <l>(X), A<0. 

Where the first (second) half-circle contour integral in the above equations is 
carried out in the lower (upper) half complex plane. The functions 4>±{ijo) are 
analytic in the upper (+) and lower (-) halves of the complex plane. Applying 
the forward Fourier transformations then gives 

+£rH = <M^)+0-H (68) 



a + {u])a-{tjj) 



where 



a + (u)a-(u) = [1 - CM]" 1 = 6 ' (69) 

2 cosh ^ 

Explicitly, 

M«) = M-) = ^r(^)" 1 7r(l-^) (TO) 

^ - - - ■ ^ 

Here T(z) is the gamma function. The common strip of analyticity for a±{oo) 
is uj G (— ivr, m). 

In the above equations, a±{oj) factorize the kernel and are analytic in the 
upper and lower halves of the complex plane, respectively. From equation (f68|) . 
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we obtain 



= + (f)-(<jj)a-{uj) - y\ (u)a-(u). (72) 



Here we have used the decomposition <p + (u)a^(uj) = $+(uj) — where 
3>+ is given by (the evaluation of <£_ is not necessary in the following) 

^-.f^."!'--^-*)). (73) 
+ v ; I 2(w + ie) w + ivr / v ; 



Owing to the analytic properties of both sides of (|72j) . they should equal some 
entire function, since both sides tend to zero as u — ► oo. Here this entire 
function turns out to be zero. 
At the Fermi point 



yi 



( ) = e (i)(B) = _ jf ^(^ = (74) 



where yi(cu) is analytic in the upper half of the complex plane. Therefore 
lim^i^oQ ujyi{uj) = holds, which yields 

e- B ' = /f-/"? ■ ■ (75) 

47TJ|| a_( — 17TJ 

We see that the Fermi boundary decreases monotonically with increasing mag- 
netic field. Subsequently, we have 

y+H = $+Ha+H = iij cff (-L J—) a-(p)a+(u). (76) 

\ a; + le w + nr / 

Similarly, y~-[oS) can be obtained from the right hand side of equation (fT2|) . 
Nevertheless, we only need y+(uj) for calculating the free energy. From equation 
CD) we then have 

ai(u)ln(l + e" £ i ^/ T )(fo 

-oo 



1 

2?r 



+ — / e-H/2 



OO 



G(w)-1 e (1)+ (u;)cL; 
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i roc fW+fuj) 

= m o) + i(i- w9 i/)--/_^^U 

= F(0,0) + ^(J 1 - M9 H)-(^) 2 2^±M. (77) 

Here finally a_(0)a+(0) = 2. The zero field free energy F(0,0) = -J|j(^(l) - 
^(2)), where ^(a) is the digamma function. It follows that the susceptibility 
behaves like x ~ 4^2 9 ] , indicating su(2) symmetry at this point as the field 
diminishes the rung interaction. This susceptibility is consistent with that of 
the isotropic Heisenberg XXX chain, where Jn = | and h = H/2 [5]. It is seen 
that at the inflection point the model can be mapped onto the su(2) isotropic 
Heisenberg chain with zero field. 



4 The quantum transfer matrix and nonlinear integral equation 
approach 

In this section we review the recently developed high-temperature expansion 
(HTE) formalism for integrable models, including the quantum transfer matrix 
(QTM) approach and the derivation of the temperature-dependant free energy 
via nonlinear integral equations with only a finite number of unknown func- 
tions. This paves the way for the application of this approach to the integrable 
spin ladder models and to the physical description of the strong coupling lad- 
der compounds. The final results can be applied without any knowledge of the 
QTM and HTE. This can be done by inserting the chemical potential terms \i{ 
of a given model into the free energy expansion (|128p . which gives instant ac- 
cess to finite temperature physical properties, as shown for a number of ladder 
compounds in section [5j 

The QTM relates a (i-dimensional quantum system at finite temperature to 
a (d+ l)-dimensional classical system on an inhomogeneous lattice via the so- 
called Trotter-Suzuki mapping [135U136tll37pi38lll39pi40] . In one-dimension, 
it involves a general relation between the transfer matrix t(u) = e lP + u,H +0(u-) 
and the Hamiltonian 7i of a quantum system, where P denotes the momentum 
operator. One also introduces an auxiliary transfer matrix t(u) which is adjoint 
to t(u) such that t(u) = e -lP + MW +°(""). In this way the partition function Zl 
of the one-dimensional system at finite temperature is mapped to that of a 
staggered two-dimensional vertex model, with Zl = tr e~^ w = limjv^oo Zl,n, 
where = l/T and Z l ,n = tr [t{-(3 /N)i{-(i /N)] N / 2 . The trace tr is carried 
out on the quantum space, as opposed to the auxiliary space. In the following 
we show explicitly that the partition function of the one-dimensional quantum 
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system is given by Z^^ = tr [£qtm(0)] l , where the QTM tQTM(^) is defined 
on the columns of the two-dimensional lattice. The QTM forms a commuting 
family, i.e., [£qtm (u) , £qtm {v)] = 0, by the virtue of the Yang-Baxter equa- 
tion. Thus the problem of describing the thermodynamics of one-dimensional 
quantum systems reduces to finding the largest eigenvalue of the QTM of a 
finite two-dimensional classical lattice model. In particular, the free energy per 
unit length is given by 

/ = - lim -3-lnZ LjA r = -^lnA max (0). 

L^oo Lp p 

In general the next leading eigenvalues of the QTM give the exponential cor- 
relation length £ of the equal time correlations. 



4.1 Preliminaries 

Our starting point is the su(r+l) Uimin-Lai-Sutherland, or permutator, model 
[9"lllUpil| , for which one can define the QTM |135pi36pi37pi38pi39pi4Uj and the 
T-system |144^145|ll46j . A description of the QTM approach for the Uimin- 
Lai-Sutherland model can be found in |141|jl42] . This involves working with 
the two-dimensional classical counterpart of the Uimin-Lai-Sutherland model, 
which is a rational limit of the Perk-Schultz model |14j . 

Let {|1), |2), . . . , \r + 1)} be an orthonormal basis of an r + 1 dimensional 
fundamental representation V of su(r + 1), and {(1|, (2|, . . . , (r + 1|} be its 
dual, with (a\b) = S ao . The i?-matrix is given by [14] 

r+l r+l 

R(v)\a) ®\b) = J2H !«') 8 \ b ') R ab\v) (78) 

a'=l b'=l 

with the matrix element (see figure [5]) 

Rab'(v) = (a' I ® (b' I R(v) \a) ® |6) = v 5 aa ,5 w + 5 a , b 5 ab/ (79) 

where v G C and a, b, a', b' £ {1, 2, . . . , r + 1}. 

Further let Vi, V 2 , . . . , V L , V ± , V 2 , . . . , V N be copies of V, i.e., Vj = Vj = V, 
where L,N/2 £ Z>i. The row-to-row transfer matrix with periodic boundary 
conditions, which acts on the quantum space V\ (8) V% <8> ■ ■ ■ <8> Vl, is defined as 



t{v) =tVyjR 2k:L {v)---R2k,2(v)R2k,i(v)), k = 1,2,..., N/2 (80) 
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Figure 5. The Boltzmann weights R^'^iv) of the vertex model. 
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Figure 6. The transfer matrix t(v)\ c t 1 '"' , " L } . The dots denote summations in definition (1821 , 

v '{Pl,---,PLI 



where R a b(v) acts on V\ ® Vi ® • • • ® Vn ® V\ ® V2 ® • • • ® Vl such that 
Rab(v)\vi) ® • • • ® |f a ) <8> • • • ® |^at) ® (8) ■ ■ ■ ® \a b ) ® • • • ® |a L ) 

r+1 r+1 

= ^ ^ |^i) ® ■ ■ ■ ® |£ a ) ® • • • ® |^iv) ® K> ® • • ■ 



5a=l»=l 



® |?76}® ■■■® \a L )Ril' v a b b (v). 



(81) 



Here Vj,Oij £ {1, 2, . . . , r + 1} and the matrix element is defined as (see figure 

ED 



^UZTJ = ^ ^ ' ' ' ® l/%> ® • • • ® W (82) 



where ctj,(3j £ {1, 2, . . . , r + 1}. 

The Hamiltonian H = TLq + TC c h of the model, which acts on V± ® V2 ® 
• • • ® Vx, consists of two parts. The first part "Hq characterizing the intrachain 
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interaction is given by (see also section [5]) 

H = J { Slnt(v)\ v=0 . (83) 

The second part is of the form 7i c \ Y = 5Zj=i^j> wnere T~£j acts non-trivially 
on the jth site of the quantum space Vj. It contains the chemical potential 
terms which include the rung interaction and magnetic field. Explicitly, 

L L 

H ch = J±J2 Si ■ ?j ~ VB9H J2( S j + T D ( 84 ) 

j=l 3=1 

for Hamiltonian (|3|) . As discussed in section 12.21 the chemical potential terms 
can be diagonalised in a suitable eigenbasis. The chemical potentials \i a are 
given by 

(ai | ® ■ ■ • ® (a L | Hj | fa) ® • • • ® | @ L ) = fi aj S aif3l ■ ■ ■ 5 aL p L (85) 

for j = 1,...,L. For example, for Hamiltonian ([3|), the eigenbasis can be 
chosen as given in equation ([5]), for which the chemical potentials are given by 

m = fj, 2 = -HBgH, ^ 3 = o, fi4 = fi B gH. (86) 

We also define an auxiliary transfer matrix with periodic boundary condi- 
tions, which acts on V\ ® V2 ® ■ ■ ■ ® Vl as 

t(v) = tT %k i (R 2k -i,L( v ) ■ ■ ■ R2k-i,2( v )R2k-i,i(v)), k = 1, 2, . . . , f . (87) 

Here R(v) is essentially a 90 degree counterclockwise rotation of R(v), whose 
components are given by =-R& 1, ^(w) (see figureE]). R a b{v) acts on the 

space V\ ® V 2 ® • • • ® V N ® V\ ® V 2 ® • • • ® V L such that 

R a b(v)Wi) (8) • • • ® |f a ) (8) • • • ® |^jv) <8> |ai) (8> • • • <8> |a&) (8> • • • ® 

r+l r+1 

= X) I] h)®---® IU® ••• 

®\v N ) ® |ai) ® • • • ® |%) ® • • • ® |a L )i?^'2>) (88) 
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Figure 7. 90 degree rotation R a , 
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Figure 8. The transfer matrix t(w)j^*' defined in (1876 . Summations are taken at the dots. 

where Uj, ctj £ {1, 2, . . . , r + 1}. The matrix element is defined as (see figure [8]) 

= <«ii ® • • • ® («£it(«)iA> ® • • • ® (89) 

where ay, /3j G {1, 2, . . . , r + 1}. 
The logarithmic derivative of (187j) also reproduces Wo- Explicitly, 

t(v) = t(0){l + H v/J {l +O(v 2 )} 

t(v) = F(0){1 + Wo u/J|| + 0(w 2 )}. (90) 

As t(0) and i(0) are, respectively, left and right translation operators, then 
t{y)t{v) = 1 +2H v/Ju +0{v 2 ). (91) 
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4.2 Quantum Transfer Matrix method 

The partition function Zl can now be written in terms of the transfer matrices 



as 



Z L = tr Vl ®-®v L e~ m = tr Vl9 ... 9Vli e- fnu e- pHo 

= Hm tr Vl8 ...^e-^» ( 1 - 

>oo y iv y 

= lim trvig-aVbe-^Ct^jv)^))" (92) 

where iV is the Trotter number, defined by utv = —J\\/3/N, with /? = 1/T the 
inverse temperature. Here we have set the Boltzmann constant to be unity 
and have used (|9ip , along with commutativity of TIq and Ti^ and the formula 
limAr^oo tr (1 — jf) N = tr e~ A for a suitable operator A. The free energy per 
site is 

/ = - lim -i- In Z L 

L^oo Lip 

= - lim lim ^\ntY Vl ^...^ VL e~ f5HA {t{u N )t(u N ))^ . (93) 
L— >oo iV— >oo pL 

In general, the eigenvalues of t{u^)t{u^) may be infinitely degenerate in the 
Trotter limit N — > oo (un — > 0) and taking the trace in (I93h is not a trivial 
operation. To overcome this difficulty, we define an inhomogeneous transfer 
matrix - the QTM - with periodic boundary conditions, 

*qtm(w) = tr^.e~ /3Wj R N j(u N + iv)R N ^ij(u N -iv)--- 

R4j(un + iv)Rz,j{uN — iv)R2,j{uN + iv)Rij(uN — iv) (94) 

which acts on V\ <8> V2 ® ■ ■ ■ <8> Vjv- The matrix element of the QTM is defined 

as 

W u )fcX} = <«i I ® • • • ® («iv i *qtm(*oia) ® • • • ® i/%> 

where ctj, (3j G {1, 2, . . . , r + 1}. Explicitly, it may be written as (see figure [9]) 

*qtm(^{^:;;; } } = 



38 



M.T. Batchelor, X.-W. Guan, N. Oelkers and Z. Tsuboi 



b N 



UN — IV 
&JV-1 L 



1 aN 



UN — IV 

bi <■ 



un + iv 
^ a 2 



a i 



Figure 9. The quantum transfer matrix tQTM (f ' bw) The dots denote summations in 
definition ((95j . The large dot stands for e _fl7i J , j 6 {1,2,..., L}. 



N/2 

{u k } k=l 



where i^tv+i = and z/& G {1,2, ... ,r + 1}. The QTM is connected to (see 
also figure fT0|) 



trv l9 ... 9 v L e ^ nch (t(uN)t(u N ))2 

N/2 

= trv l8 ...®vte" ^^^ n^^-i^^J^fe^^iv) • • • R2k,2{UN)R2k,l{UN)) 

k=l 



X (R 2 k-1,L(UN) ■ ■ ■ R2k-l,2{uN)R2k-l,l{uN))} 
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► 4 

u + iv 


\ i 

u + iv 
~\ 


> i 

u + iv 

-\ 

—i ■ 


> 

u + iv 


— iv 

r 


u — iv 

r 


u — iv 

r 


u — iv 

r 















I 





u + iv 


u + iv 
^ . 


u + iv 


_u + iv 


— 

r 


u — iv 

r 


u — iv 

r 


u — iv 

r 















Vi V 2 Vl-i V l 

Figure 10. The QTM lattice structure. Summations are taken at dotted points. Large dots stand 
for e - * 9 ^ , j e {1, 2, . . . , L}. The values u = njv and v = correspond to equation |j96j . Horizontal 
lines with u + iv correspond to t(v) while horizontal lines with u — iv correspond to t(v). The 
vertical lines correspond to iQTM C^)- 



N/2 

tlV l8 ... (8 v I1 e~ ^UPHj JJtly 2ki(S) y 2k {(R2k,L( u N)R2k~l,L( u N)) ' ■ ■ 
k=l 

*{R2k,2{uN)R2k-l,2{uN)){R2k^{uN)R2k-l,l{uN))} 
V NJ2 

tr v 1 ®---®v N Tl tl v 3 e~ m3 n(R2kj(uN)R2k~ij(uN)) 

j=l k=l 

tr ^.WW0)) L - (96) 
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Here the ordered product is denned by 



L 

l[A(j) = A(L)...A(2)A(l) 

3=1 

for any indexed operator {A(j)}. 

Let Ai, A2, . . . , A( r+1 )N be all eigenvalues of £qtm(0) such that Ai > A2 > 
. . . > A( r+1 )« . Then the free energy per site (|93|) is 

/ = -lim lim J_InAf{l+f^V+f^V + 



A( r+1 y 



L " 



It is assumed that there is a finite gap between the largest and the second 
largest eigenvalues, Ai and A2, of £qtm (0) , at least when A" is a finite number. 
This gap is expected to still exist in the Trotter limit if T > 0. Moreover, 
Ai, A2, • • • , A( r+1 )N do not depend on L. Thus if we exchange [135|,ll36(.ll38j 
the order of the two limits L — > 00 and N — > 00 in (197h . the free energy per 
site can be expressed in terms of only the largest eigenvalue Ai of the QTM 
flSD at v = 0. Namely, 

f = -T lim log Ai (98) 

N—>oo 

a remarkably simple result. 

Now the eigenvalue [EElUmUBI] T^\v) of the QTM ®ty derived by the 
Bethe Ansatz is 

r+l 

T^(v) = ^2z(d;v). (99) 
d=i 

The functions {z(d;v)} are defined as 

z(a; v) = Mv) — ; —, — ( 10 °) 

Qa-i{v - 51(0 - l))Qa{v - 51a) 
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for a G {1, 2, . . . , r + 1}, with 

Q»=n^-4 a) ) 

fc=i 



where M a G Z>o and QoO ) = Qr+i(v) = 1- The vacuum parts ijj a (v) appear- 
ing in (jlOOp are given by 



4W=e""" / ^ + (" + iV)^("- i ^i) ( 101 ) 

for a G {1, 2, . . . ,r + 1} and 0±(u) = (v ± iujv)^ 2 . The complex numbers w|, 
satisfy the Bethe equation 

wrf + _ ^-i(4 a) + ii)Qa(4 a) - i)Q a+ i(4 a) + & (102) 



2 1J 



for k G {1, 2, ... , M a } and a G {1, 2, ... , r}. 

All of the above is for general r. Specializing to the su{4) (r = 3) case of 
relevance to the integrable ladder model ([3j) , with chemical potential terms 
([86]) . we have 



-c 



+ Q 1 (4 1) +i)Q 2 (4 1) -ii) 
x = c ^ Qi(4 2) + gQgrf - i)Q 3 (4 2) + gi) 

<M4 3) + fi) _ ^ ^(4 3) + li)Qs(4 3) - O 



(103) 



-c 



<M4 3) + fi) " Q 2 (4 3) -|i)Q3(4 3) +i)' 

In the first line, k G {1, 2, . . . , Mi}, in the second, G {1, 2, . . . , M2}, and in 
the third k G {1, 2, . . . , M 3 }. 



4.3 T-system 

The next step is to introduce an auxiliary function Tm\v), whose origin is 
related to a fusion hierarchy |177| of the QTM (cf. the Bazhanov-Reshetikhin 
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formula in [178] ). 



a m 

T^(v)= ]T [J n z(d jtk ; v-±i(a-m-2j + 2k)) (104) 
K, fc } j=ik=i 



where the summation is taken over dj t k £ {1, 2, . . . , r + 1} such that dj t k -< 
dj+i,k an d djjt ^ dj t k+i (1 -< 2 -<■••-< r + 1), with m G Z>i and a £ 
{1,2, ... , r}. This general function Tm\v) contains the eigenvalue T^\v) in 
(j99|) as a special case, namely (a, m) = (1, 1). The poles of Tm (v) are spurious 
due to the Bethe equations (I102j) . The relation (|104|) can be interpreted as 
Young tableaux with spectral parameter v if we set z(a; v) 
For example, for the sit(4) case (r = 3), we have 



a: v 



U) 



hv 



+ 



2;v + 3;v 



1 


f — 


* 


2 


u + 


r 



+ 



4: f 



(105) 



3;v + p 



+ 



l:v 



4;v + ±i 



+ 



2; 


u — 




3; 


u + 


? 



+ 



2;v-M 
TT 



4;^ + ±i 



3; v 



4;-u + ±i 



(106) 



l;u + p 


l;v- p 


+ 


l;^ + p 


2 ;u -fi 



+|l;u + p 



4;u 



2;v + p 



4;u- ii 



2;t> + |j 



2;t; - ^i 



3: v + li 3: v 



l;v + p3;t> — p 



2;u + p 



3; w 



3; v +"p 4; u — ~i 



4;u + p 



(107) 



Here we have used the notation 



z(dix,v- h)z(d 2> i;v + p) 



^2.15^ + ni 



z(di 5 i; v + p)z(di >2 ; w - p) = \diy,v + p|d 1;2 ;?J - p r 



In this sense, relation (|104p represents the character of the m-th symmetric 
and the a-th anti-symmetric tensor representation of su(r + 1) with spec- 
tral parameter. For a G {1,2, ... ,r} and m E Z>i, it may be normalised as 
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(«) = T^{v)/N^{v), where 

- ii(a + m))(j) + (v + ii(a + to)) 



0-(i> — |i(a — m))(f)+(v + ii(a — to)) 



x ] J - ±i(a - to - 2j + 2A;))0+(u - ±i(a - to - 2j + 2jfe)). (108) 

j=l k=l 

One can show that Tm\v) satisfies the functional relation 

f + |i)f - ii) = r^iM^ («) + fir 1 ^)^ 1 ^) (109) 



which is called the T-system |144|I145|I146| . Here a £ {1,2, ... ,r} and m G Z>i, 
with 



T (a) (i;) = l, for aGZ>i 
~ {0) 0_(r; + iim)^ + (w-^im) 



4>~(v — | im)(j)+(v + ^ im) 

-m(^ 1 +^ 2 H hAV+l) 

= e r for to G Z>i. 



(110) 



It can be verified, for example, that (|105p - (|107p satisfy the T-system (|109D 

4 4 

T x (1) (t; + lijT^Cv - ^i) = E E 6Z±H^EB 



a=l 6=1 



a<£> 



a;v + \\ 


b; v — \\ 




a; v + \i 


b; v — \i 



a>b 



a<b 
,(1), 



a>b 



b; v — ^i 



a; t> + ^i 



T^{v)+T[ 2 \v). 
Then taking note of the relations 



(111) 



M[ X \v + \i)N[ X \v - \i) = N^\v) = jVi 2 \v)/fi 0) (v) (112) 
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gives the desired result 

+ Ji)if >(„ - ±i) = f£\ v ) + f[°\v)f[ 2 \v). (113) 

4.4 Nonlinear integral equations 

The thermodynamics of the spin-^ Heisenberg XXZ chain has recently been 
described by a nonlinear integral equation (NLIE) |179| . which was rederived 
[180 J from the T-system j!46j of the QTM. It was also derived from a fugacity 
expansion [181j . There are similar NLIEs for other models |130|.I1311I147|,I148| 
149J. Most importantly, the number of unknown functions for the NLIE is 
only one, corresponding to the rank of su(2). In this subsection, we review the 
derivation of the su(r+l) NLIE (|123p from the T-system (|109p . in which only r 
unknown functions appear. This paves the way for the exact high temperature 
expansion to follow in section 14.51 

Numerical analysis for finite values of N, um and r reveals the expectation 
that the largest eigenvalue of the QTM (|94p at v = is given in terms of 
a one-string solution (for each colour a = 1,2, ... ,r) of the Bethe equations 
(|102p in the sector N/2 = M\ = M2 = . . . = M r . From now on, we consider 

only this one-string solution. Thus T± (0) should give the largest eigenvalue 
Ax- Numerical evidence for this conclusion is given in Appendix C. We thus 
state the following conjecture. 

Proposition 4.1 For small un (\un\ <C I), a G {1,2, . . . ,r}, all zeros {z^f 1 } 
of T^ a \v) are located outside the physical strip Imv E [— iji]. In particular, 
they are located near the lines Imv = ±^(1 + a), at least for the case fi\ = 
f>2 = ■■■ = Mr+1 = 0. 

This conjecture is expected to be valid even in the Trotter limit N — > 00. 
Now Tm\v) has poles only at ±/3m\ where [3$ = \i{m + a) + iujv, whose 
order is at most N/2. Moreover, 

= lim f[ a \v)= V e -^e"^ ... e~^r (114) 

\v\ — >oo — J 

l<ii<i 2 <-..<i a <r+l 

is a finite number. This quantity is related to a solution of a functional relation 
called the Q-system [IB21IIB3] . for which (Q^) 2 = Q^Q^+i + Qt'^ Qt +l \ 
where m £ Z>i and a £ {1, 2, ... , r}. 
In the su(4) case, we have 



Ml M2 >'3 _£4 
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(2) Ml +M2 Ml +M3 Ml +M4 M2+M3 M2+M4 M3 + M4 

(3") M1+M2 + M3 Ml +M2 + M4 M1+M3 + M4 M2+M3+M4 

= e T +e T +e T +e T 
Q (4) =e - - + ^ +> ' 4 , (115) 

It follows that we must put 

JV/2 ( Aa) -r(a) } 

Tl a) (v) = Q<°> + V { 3 —— + \ (116) 

where the coefficients b^ , b^ £ C are given by the contour integrals 



b ? - w v+ & y - (n7) 

c 



The contour is a counterclockwise closed loop around (3[ which does not 

encircle — f3[ a \ Similarly the contour C is a counterclockwise closed loop 

around — p[ a ^ not encircling p[ a \ 
Using the T-system (I109p . the integrals (|117h can be modified to read 

b (a )= f dv f f^\v-li) 



^\v-\i)T^{v-\i) 



f[ a \v-i) 



+ M ; - ' ^ K« + (118) 



%w 2m { T i a \v + i) 
fi a) (v + i) 

Substitution of (|118p into (|116p and performing the summation over j then 
gives 

f[ a \v) = Q^ 
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+ 



/ \ N/2 



cw 27ri v-y- p{ a) I T[ a) (y + p[ a) - i) 

ft'Hy + # o) - W[ a+l \y + (3{ a) - ii) 



+- 



+ 



N/2 

dy 1 '^) {T!f\y-^ ) + ^ 



C<"> 27Ti y _ y + I T (a) {y _ £(a) + j 



T w (y _^) + i) J 

for a G {1,2, ...,r} and m £ Z>i. Here the contour C( a ) (resp. C ) is a 
counterclockwise closed loop around the origin which does not encircle — 2(3^ 

(resp. 2p[ a) ). There are cancellations of singularities in the denominator and 
the numerator in each term in the first curly bracket {• • •} in (1119p . The poles 
in the first term in the first bracket are -2(3{ a) and z^ - (3{ a) + i; the poles 

in the second term in the first bracket are 0, -2f3[ a) and z[ a) - f3[ a) + i. Thus 
the contribution to the contour integral from the first term in the first bracket 
vanishes if the contour does not encircle the pole at z^ — p[ a ^ + i (cf. 
Conjecture 14. ip . As for the second curly bracket in (j!19|) . the poles in the first 
term are 2/3^ and + p[ a ^ — i; the poles in the second term are 0, 2^ 
and + p[ a ^ — i. Thus the contribution to the contour integral from the first 

term in the second bracket vanishes if the contour C° does not encircle the 
pole at zf + /3( a) - i. 
It therefore follows that (|119p reduces to 

Tl a) (v) = Q<f ) 



+ 



sN/2 

dy 1 - Ujfe) Tt 1] (y + /3[ Q) ~ ^)Tj a+1) (y + ffl - £i) 
CW 27ri v-y- p[ a) Tl a) {y + $ a) - i) 



/ xN/2 

+ I dy 1 - ft 1] (y - P[ a) + ^)fj a+1) (y - Pi a) + ji) fl2Q1 

for a € {1,2, ... ,r}. Here the contour (resp. ) is a counterclockwise 
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closed loop around the origin not encircling z[ — (3[ + i, —2/3^ (resp. + 
(3{ a) - i, 2f3~{ a) ). The poles at y = in f[ a ~ l \y ± /3| a) =p |i) are cancelled by 
the zeros at y = in (^p)^ 2 - Thus ([T2D]) is further simplified to 

f[ a \v) = Q { ? ) 

dy Tt x \ y + ^ - \0[ a+1 \y + p{ a) - ±i) 



+ 



oW 2vri (« - y - ^T?' (y + /3( a) - i) 



27Ti („ _ y + pWft]*) { y _ ft) + J) 

for a € {1,2, . . . ,r}. 

This is the NLIE for finite Trotter number. The next important step is to 
take the Trotter limit N — > oo. Define T^ a \v) = YmiN->ooT^ (v), then 

Finally, we have arrived at a system of NLIE which contains only a finite 
number of unknown functions, {T^ a \v)}i<. a <. r , namely 



Ti a \v) = + I 
Jc 



4 dy ^% + /?j a) - foTj^fr + /f } - |i) 

CC«)27Ti (v _ y _ (3 (a) )T (a) {y + f} (a)_ i) 



Jc'»> 2vri (u _ y + ft))7< a \y - (3{ a) + i) 

for a £ {1, 2, . . . , r}. Here T^ r+1 \v) = exp(— {yt\ + • • • + /i r+ i)/T). The contour 

(resp. C (a) ) is a counterclockwise closed loop around the origin such that 

y / v — p[ a ^ (resp. y ^ v + (3^) and it does not encircle z[ a ^ — ft^ + i, —2(3^ 

(resp. Zi + /?| a ' > — i, 2/3^ ), where 2^ = limyv->oo -^j' 1 ''. 
In particular, for the su(4) case, the NLIEs are 

Jy W (2) t 
ICW*™ (v - y - (y) 



/en) 2vri fw-w-iVrWfV) 
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+ 



dy ^(wkr) 1 } 2 ^-^ 



+ 



C (2) 27T1 ( V _y_* i ffi*)(y + |i) 

U« 2 > 2vri ( v _ 2/+ | i) ^(2) (2/ _ i i} 



Ti (3) w=Q (3) +e - » 1 ^ r ^ /• 

/CO) 2tti (y-y- 2i)Tf'(y + i) 



Jc^ 27ri( u _ y + 2i)T 1 (3) (y-i) 

where are given in (|115p . 
The free energy per site follows from (|123p and the relation 

/ = Jr. -TlnT/^O). (125) 



4.5 High temperature expansion 

We now turn to the calculation of the high temperature expansion of the free 
energy (|125|) from the NLIE (|123|) . In the su(2) case (the isotropic Heisenberg 
chain), the high temperature expansion of the free energy was calculated from 
the NLIE |179| up to order 100 |129j . The first step is to assume the expansion 

ri a) (v) = exp (f>8°(«) (^f) n )< «e{l,2,...,r} (126) 

for large T. In contrast to the su(2) case |129j . poles in the coefficients {bffl(v)} 
have to be taken into account. This requires the further assumption, 

n-l c ^ v 2j 

^(«)=E ~TTlV> (127) 

where the coefficients c^- do not depend onu, In general, these coefficients 
are rational functions of Q± given in (|114|) and thus they are T-dependent. 



Integrable quantum spin ladders: comparison between theory and experiment 



19 



The coefficients bn (v) are obtained by substituting (|126j) into (|123|) . with 
= log Q< a >. 

For the sti(4) case relevant to the integrable ladder model the first few terms 
of the free energy are given by [131] 

- l -m H) = In 0« + c« (£) + cg> (£) 2 + c« (£) 3 + • • • (128) 

The coefficients c£ q are given up to order n = 5 in Appendix A. The key 

point for the ladder model is the quantity given in (1115j) which involves 
the chemical potentials. 

The advantage of this approach is that the thermal and magnetic properties 
can be evaluated directly from the exact free energy expression (I128j) using 
the standard relations 

df(T,H). d 2 f(T 1 H) i „ d 2 /(r,.H~) , , . 

M = — - t y = v — - t C = -T — v — - r (129) 

which define the magnetization, the magnetic susceptibility and the magnetic 
specific heat. 



5 Thermodynamics of the spin-- ladder compounds 

Although many two-leg ladders have been synthesised and extensively studied, 
the critical points have been measured only for a few strong coupling com- 
pounds with accessible gaps. Theoretical results for the two-leg spin-^ ladder 
model have been mainly focussed on the thermal and magnetic properties of 
the standard Heisenberg ladder. According to a celebrated result from pertur- 
bation theory [63 P 64 1 I20]. the first-order terms for the zero temperature gap A 
and the critical field H C 2 are given by A = J± — Jm and [iBgHci = J± + 2Jy, 
which are seen to be in good agreement with the experimental results for some 
strong coupling compounds [20J. In addition, the susceptibility as a function 
of temperature, % ^ e~^^ T /VT, derived from the quantum transfer matrix 
algorithm [66] is generally used to fit the susceptibility of the two-leg ladder 
compounds. 

In comparison with perturbation theory or numerical and analytic meth- 
ods [3" Hll62l[163l l2"T lll84j applied to ladder-like systems, the integrable ladder 
model ([3]) [92 ^11 16^1 17t l8| [T52"] has a distinct advantage for studying phase di- 
agrams and quantities such as the high field magnetization, specific heat and 
susceptibility, as they follow directly from the exact free energy. Moreover, the 
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Coixip ounds 


B5i2aT 


Cu(Hp)Cl 


BPCB 


BIP-BNO 


CuCON 


HSL J_l 


13. OK 


13.1 K 


13.3 K 


72 K 


509 K 


HSL J|| 


1.15K 


2.62 K 


3.8K 


17K 


44 K 


ISL J± 


13.3 K 


14.5 K 


15.4K 


83 K 


520 K 


ISL J,, 


0.29 K 


1.0K 


1.2K 


7.0 K 


10K 


HSL H c \ 


8.3 T 


7.5 T 


6.6 T 


38.8 T 


323 T 


TBA H cl 


8.6 T 


7.4 T 


7.4 T 


41 T 


334 T 


HSL H c 2 


10.4T 


13.0 T 


14.6 T 




415 T 


TBA H c2 


10.02 T 


13.1 T 


14.1 T 


83 T 


389 T 



Table 1. The coupling strengths, J ± and J\\ , and critical field values, H c \ and H c 2 obtained from previous 
fits with the Heiscnberg spin ladder (HSL) and for the intcgrablc spin ladder (ISL) using the TBA approach 
for some strong coupling ladder compounds (see text). 



ground state properties at zero temperature are now understood by means of 
the TBA |1161I117] . As demonstrated in section [31 the critical fields, H c \ = 
(J_i_ — &J\\)/ hb9 and H C 2 = (Jj_ +&J\\)/ iib9, also give a good fit for the strong 
coupling compounds. The coupling constants J± and Jn can be fixed by fitting 
experimental results for the susceptibility and the magnetization. In general, 
these coupling constants cannot be fixed from a good fit just to the susceptibil- 
ity, as it determines only the amplitude of the gap. The way to determine the 
coupling constants has been discussed at the begining of section [331 In general 
the values suggested from only the critical fields are not the best for fitting. 
The coupling constants should provide a good overall fit for each of thermo- 
dynamic properties. In the following we apply the integrable model approach 
to the compounds (5IAP)2CuBr4 • 2H2O [42] . bis 5-iodo-2-aminopyridinium 
tetrabromocuprate (II) dihydrate (abbreviated B5i2aT), Cu2(C5Hi2N2)2Ci4 
(abbreviated CuHpCl) [55153155] . (C 5 Hi 2 N ) 2 CuBr 4 (abbreviated BPCB) [47], 
BIP-BNO [IMKIHO] and [Cu2(C 2 02)(CioH 8 N2)2)](N0 3 )2 (abbreviated Cu- 
CON) [THT] . 

Values for the critical fields, obtained from the integrable spin ladder (ISL) 
via TBA at zero temperature |116]lll7] . are compared with the experimental 
results for the Heisenberg spin ladder model (HSL) at low temperature in table 
1. The small discrepancy between the TBA critical fields at zero temperature 
and experimental critical fields at low temperature is a finite temperature 
effect. As the temperature increases the energy gap decreases due to triplet 
excitations, as the energy gap is a measure of the energy cost of the first 
triplet excitation at zero temperature. The experimental gap is determined 
by the first critical field H c \ at low temperatures. It will be seen below that 
the thermal and magnetic properties derived from the HTE are in excellent 
agreement with the experimental results. 
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5.1 Note on conversion constants 

Before turning to the study of the thermodynamics of these compounds via 
the HTE, a remark on conversion constants is in order. Overall conversion con- 
stants for the susceptibility, the specific heat and the magnetization are needed 
due to the fact that the T-system (jlQ9j) is specified up to normalization. There- 
fore, in order to match the experimental data for the unnormalized magnetic 
properties, a scaling factor needs to be fixed for the magnetic properties, such 
as the magnetization, susceptibility, specific heat and magnetic entropy. In 
the following, we will see that the conversion constant xhte ~ 0.40615^exp 
for the susceptibility of the ladder compounds is half that of the spin-1 chain 
compounds |153^ll54j due to there being 2L spins in the ladder model. The 
specific heat conversion constant is material dependent due to the choice of 
experimental units. However, for normalized magnetic properties, such as the 
magnetization (M/M s ), the conversion constant plays no role in fitting the 
experimental data. 

From the derivation of the HTE it can be seen that the free energy expansion 



Physically this amounts to the fact that the units of energy, temperature, ex- 
ternal field and coupling strength are not unified in the theory, i.e., there is no 
analogue of a 'characteristic length'. Mathematically, it can be traced back to 
the unit of temperature and the multiplicativity of the i?-matrix of the model. 
The scaling factors for different magnetic quantities such as susceptibility, spe- 
cific heat and magnetic entropy can be uniquely fixed by their corresponding 
experimental curves. They depend on their experimental units and are given 
in the following figure captions. 

In the HTE calculations for physical properties this factor has previously 
been used to obtain a best fit for experimental data via vertical scaling. This 
procedure is open to critique, as it would offer a new free parameter, enabling 
a better fit than from the pure HTE theory. The conversion factor can in fact 
be calculated, leading to predictions of physical properties in physical units 
(SI) instead of arbitrary HTE units. One remaining problem is the widespread 
use of non-standard non-SI units in the experimental measurements. These 
measurements sometimes even depend on compound-specifics such as molec- 
ular mass. The problem of converting SI units to experimental units, while 




(130) 



is only fixed up to a normalisation constant, i.e. 



/phys. units — Tconv. ■f H TE(T,H). 



(131) 
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often not trivial, is in no way related to the HTE method and will thus not 
be further discussed here. HTE predictions are made in SI units. 

In eqns (|130p and (|13ip the unknown constant ^ C onv. is a pref actor, i.e., the 
same in all orders of the expansion. The expansion is generated for interacting 
spins from the Bethe Ansatz solution. For zero order, the non-interacting case 
J = 0, the expansion reduces to the trivial (non Bethe Ansatz) case of L 
independent spins (here, e.g., a two site spin-^ system), for which 



The second equality holds for the special point J = 0, but since 7 con v. is a 
constant it can be extended in eqns (|130p and (|131j) to all J > 0. 

Once the free energy is thus known in SI units, it is straightforward to obtain 
all other properties derived from it in SI units via 



The important point here is that the HTE method does make predictions in 
physical units, even if the original T-system is only fixed up to normalisation. 
Thus all conversion factors can in theory be calculated using the additional 
information from the non-interacting spin system or directly from fitting ex- 
perimental curves. 

In practice the HTE scheme is unable to access thermal and magnetic prop- 
erties for temperature scales less than Jy. At high temperatures the first few 
terms in the free energy expansion are sufficient to determine the thermody- 
namic properties. The HTE approach works best in practice for small ratios of 
Jy /T, as only the lowest few terms of the series are needed to get a good agree- 
ment with the experimental results. Consequently for high temperatures and 
materials with a small coupling constant Ju even the zeroth order HTE suffices 
to capture the main physical properties. This lowest order corresponds to the 
approximation where Ju = 0, sometimes called the rung-dimer model j!85j . In 
this case the ground state wavefunction is a simple product of rung singlets. 
To this order the partition function of the ladder model consists of only four 
terms (see the first term of the Q-system for comparison) and can be evaluated 
directly without the use of the integrable theory. However, not all interesting 
physical properties can be captured by this simple dimer model, which is the 
zeroth order of the HTE series. Especially as much of the physics of interest lies 
in the lower temperature regime, with the parameter Ju material dependent 
and not necessarily small enough for the dimer approximation. 

Figure [11] shows a comparison of the magnetization curve at two different 
temperatures for the ladder compound (CsH^N^CuBr,}. It can be seen that 




(132) 
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phys. units. 



(133) 



Integrable quantum spin ladders: comparison between theory and experiment 



53 




Magnetic field H(Tesla) Magnetic field H(Tesla) 

Figure 11. Comparison between the dimer approximation (Oth order HTE) and the 5th order 
HTE result for the magnetization of the compound (CsHi2N)2CuBr4. The dimer approximation 
does not distinguish the two critical field values, for which the experimental difference is 

H cl - H cl (see text). 

the zeroth order approximation of the pure dimer model captures the quali- 
tative shape of the curve, whereas the HTE series result (discussed further in 
figure [20)) is in quantitative agreement with the experimental measurement. 
The dimer model is obviously a bad approximation for this compound, as the 
two critical fields, namely H c \ = 7T and H C 2 = 14.1 T are far apart. The 
difference between these values is proportional to the coupling Jm , which the 
dimer model neglects. 

We proceed now to the comparison with the experimental data for a number 
of compounds. 

5.2 (5IAP) 2 CuBr 4 • 2H 2 

5.2.1 Susceptibility. Measurement of the susceptibility and magnetization 
of the compound B5i2aT [42] suggests that it lies in the strong coupling regime, 
with coupling constants Jj_ = 13.0 K, Jm = 1.15 K for the standard Heisenberg 
ladder. From the application of the HTE for the Hamiltonian ([3]) we find that 
the coupling constants J± = 13.3 K and Jm = 0.2875 K give excellent fits to 
both the susceptibility and the magnetization. The temperature dependence 
of the susceptibility curves is given in figure [T2l where the solid line denotes 
the susceptibility derived from the free energy (|128p with up to fifth order 
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Temperature T(K) 

Figure 12. Susceptibility versus temperature for the compound B5i2aT. Circles denote the 
experimental data 1421 and the solid curve is the susceptibility evaluated directly from the HTE at 

H = 1 T. A parameter fit suggests the coupling constants J± = 13.3 K and Jy = 0.2875 K with 
g = 2.1 and fis = 0.672 K/T. The inset shows the same fit in the susceptibility at low temperature. 
The conversion constant is xhte ~ 0.40615xexp- 

of HTE. An overall excellent fit with the experimental susceptibility curve 
is observed. The typical rounded peak in the zero magnetic field susceptibil- 
ity curve, characteristic of a low dimensional antiferromagnet, is seen around 
8.1 K. The susceptibility exponentially decreases with further decrease of the 
temperature toward zero. The TBA analysis suggests an energy gap A ~ 8.6 T. 

The low temperature behaviour of the magnetic susceptibility, shown in the 
inset of figure [T2l is also in excellent agreement with the experimental data. 
We observed that the magnetic field raises the susceptibility due to the lower 
triplet component being energetically favoured by the magnetic field at low 
temperature. If the magnetic field is much smaller than the temperature, there 
is no significant change to the susceptibility. 

5.2.2 High field magnetization. The behaviour of the field dependent mag- 
netization curve can lead to the prediction of the low temperature phase dia- 
gram as well as the magnetization plateaux. It follows that the magnetization 
is a particular interesting quantity for studying critical behaviour. The high 
field magnetization curves evaluated from the HTE at different temperatures 
are shown in figure[13]in comparison with the experimental curves, with which 
there is also excellent agreement. At very low temperature, provided the mag- 
netic field is less than the critical field H c i, the rung singlet forms a singlet 
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Magnetic field H(Tesla) 

Figure 13. Magnetization versus magnetic field indicating the high field quantum phase diagram 
for the compound B5i2aT with the same constants as in Fig |12l The discrepancy in magnetization 
curves at T = 0.4 K and T = 1.59 K is due to paramagnetic impurities, which result in nonzero 
magnetization in the singlet ground state. This impurity effect becomes negligible at higher 
temperatures. The inset shows the one point correlation function vs magnetic field — at low 
temperature the singlet ground state can be clearly identified. 



ground state. The antiferromagnetic correlation length is finite while the triplet 
state is gapful. At finite temperature the triplet excitations become involved 
so that the gap decreases with increasing temperature. This behaviour is ob- 
served in the high field magnetization curves at T = 1.59 K and T = 4.35 K 
presented in figure fl~3l The gap closes at the critical field H c \ = A/nsg- If 
the magnetic field is above the critical point H&, the lower component of the 
triplet becomes involved in the ground state. At zero temperature, one can 
show rigorously that the other two (higher) components of the triplet never 
become involved in the ground state. Therefore, in this vicinity, the strongly 
coupled two-leg ladder can be mapped to the XX Z Heisenberg chain with an 
effective magnetic field term, as done in section [3] and in [644167] . The mag- 
netization increases almost linearly with the field towards the critical point 
H C 2, at which the ground state becomes fully polarized. In addition, the in- 
flection point is indicated at Hip ~ J±/hb9 where the magnetization moment 
is equal to 0.5. At T = 0.4 K, the HTE magnetization curve indicates that 
H c i ~ 8.6 T and H C 2 ~ 10.02 T. These values are in excellent agreement with 
the experimental estimates of H c \ = 8.3 T and H C 2 = 10.4 T. The experimen- 
tal magnetization in the singlet ground state at low temperature appears to be 
nonzero. This is caused by paramagnetic impurities due to the synthesization 
process. 

The inflection point is clearly visible in the experimental magnetization 
curves for the compound B5i2aT [32]. It follows that for the strong coupling 
ladder compounds at zero temperature, the one point correlation function 
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(Sj - Tj) = — | lies in a gapped singlet groundstate, which indicates an ordered 
dimer phase, while (Sj ■ Tj) = \ in the fully-polarized ferromagnetic phase. 
However, in a Luttinger liquid phase (Sj-Tj) = — j J rS z and the magnetic field 
increases the one point correlation function. At low temperatures T < Jn the 

one point correlation function is given by (Sj • Tj) = \ + (37^/(^1 ^)) T ' 
field-induced quantum phase transition can be clearly seen from the one point 
correlation function curves at different temperatures (see the inset of figure 
[T3j) : in the region H < H c i, the correlation functions (Sj ■ Tj) = — | indicate 
a gapped singlet groundstate; in the region H > the correlation func- 
tions (Sj • Tj) = \ correspond to a fully polarized ferromagnetic state; while 
for H c \ < H < H C 2, the one point correlation function — | < (Sj ■ Tj) < ^ 
indicates a magnetic Luttinger phase. 



5.3 Cu 2 (C 5 H 12 N 2 ) 2 Cl4 

5.3.1 Magnetic susceptibility. The ladder compound Cu(Hp)Cl has been 
extensively studied, both experimentally and theoretically, including several 
numerical-based calculations [33U35UB51 liEMg31 M 

152J. Nevertheless, the thermodynamic quantities of the model, such as the 
high field magnetization, the full temperature susceptibility, the magnetic spe- 
cific heat, the entropy and the correlation length, are worth close examination 
via the TBA and HTE methods. Here we show that the integrable model pro- 
vides an efficient way of investigating the thermodynamics of this compound. 
The solid line in figure [TH shows the zero field magnetic susceptibility curve 
obtained by HTE with up to fifth order terms. The rounded hump is observed 
around T = 8.0 K. A full fit with the experimental curve for the susceptibil- 
ity of Hamiltonian ([3]) suggests the values J± = 14.5 K and Jm = 1.0 K with 
g = 2.1, which lies in the same phase as the conventional Heisenberg ladder 
model (64]. The TBA gap A ps 10.5 K agrees well with the experimental gap 

A = 10.9 k dHHSlES] . 

At this point it is appropriate to remark on the convergence of the HTE 
series. The HTE in equation (|126p is not a series in the usual Taylor sense. 
Here the expansion parameter is J\\/T. For a given compound Jm is fixed, 
with J\\/T a small parameter at high temperatures. It is important to note 

(k) 

that the expansion coefficients c\- depend on the temperature through the Q- 
system, which in turn depends on the chemical potentials, see equations (jl 14|) 
and (|115p . For this more sophisticated type of series expansion it suffices to 
consider only the first few terms in the physically interesting regime. When 
the expansion parameter is no longer small, adding higher order terms does 
not improve the fit significantly. In the following numerical example we show 
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Figure 14. Susceptibility versus temperature for the compound Cu(Hp)Cl. Circles represent the 
experimental data 1591 . The solid curve is the susceptibility evaluated directly from the HTE with 
/i B = 0.672 K/T, Jj_ = 14.5 K, J» = 1.0 K and g = 2.1. The conversion constant is 

XHTE ~ 0.40615xexp- 



this behavior for a generic case, namely the particular compound Cu(Hp)Cl 
under consideration (figure PH]). For the susceptibility x(T) the relative error 
in using the HTE up to third order and up to fifth order compared to using up 
to ninth order terms (as an example of higher order terms) for some sample 
temperatures is given in the following table. 



T 


2K 


5K 


10 K 


50 K 




error Xi-3 
error xi-5 


rel. to xi-9 
rel. to xi-9 


7.33% 
0.56% 


0.122% 
0.00147% 


0.0003288% 
3.8 • 10" 6 % 


< 9.2 

< 5.2 


10"''% 
10" 10 % 



In this article for comparison with experimental data we have used the HTE 
expansion with up to 5th order terms. Technically it is not harder to include 
higher order HTE terms in the analysis. The computation of the coefficients 
in Appendix |A] has to be done only once, as the expansion formula gives the 
coefficients in terms of the abstract Q-system. For each compound or new fit 
only the numerical values have to be fed into the chemical potentials, allowing 
an immediate calculation of all physical properties up to the order of the initial 
free energy expansion. 

In principle, the range of validity of the HTE can be extended to low tem- 
peratures via Pade series analysis. This uses the fact that the free energy 
expression is a smooth function of its arguments. 



5.3.2 Magnetic specific heat and entropy. The solid and dashed lines in 
figure [15] denote the specific heat curves at H = T and 4 T obtained from the 
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Figure 15. Specific heat versus temperature at different magnetic field values for the compound 
Cu(Hp)Cl with the same coupling constants J±, Ju and g as in figure H"4l The experimental 
data |59j is given by circles for H = OT and triangles for H = 4T. The solid and dashed curves are 
evaluated from the HTE for these same magnetic field values. The conversion constant is 
Chte ~ 4.515C*exp (J/mol-Cu-K). 

HTE with the same coupling constants as in figure [TH They agree well with 
the experimental data [59]. A rounded peak indicating short range ordering 
is observed around T = 4.5 K. At temperatures below T = 4.5 K there is an 
exponential decay which is believed to be due to an ordered phase. A similar 
fit for the specific heat at H = 2T and 6T is observed in figure [TBI Another 
peak is found around T = 1 K for the magnetic field H = 6T. This is mainly 
because the HTE is not convergent for this model if the temperature is less 
than 1 K. At these values the temperature scale is comparable to the coupling 
J|| = 1.0 K. Adding further HTE terms is not sufficient to overcome such 
kinds of divergency. The entropy has been calculated directly from the free 
energy (|128p with up to fifth order terms. Figure [T71 shows agreement with the 
experimental entropy curves [59] . A discrepancy with the fitting is observed 
for high magnetic field, which is probably caused by paramagnetic impurity 
effects. 



5.3.3 High field magnetization. From table 1, we see that the parameter 
fitting values J± = 14.5 K and Jn = 1.0K lead to the critical field values 
H c \ = 7.4 T and = 13.1 T, which are in good agreement with the exper- 
imental values H c \ = 7.5 T and H C 2 = 13.0 T. The experimental high field 
magnetization has been measured by different groups |64|,l65j . The zero tem- 
perature TBA magnetization curve (red dashed line) in figure [TH] indicates 
the critical field values in the proximity of the lowest temperature data at 
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Figure 17. Entropy vs temperature comparison between theory and experimental data 59_ for the 
compound Cu(Hp)Cl at the external magnetic field values H = OT, iJ = 4T and H = 6T. The 
coupling constants are the same as in the previous figures for Cu(Hp)Cl. The solid, dashed and 
dot-dashed curves evaluated from the HTE for these same magnetic field values. A discrepancy 
with the experimental entropy data for high external magnetic field is observed, probably due to 
paramagnetic impurities. The conversion constant is Shte ~ 4.1535'exp (J/mol-Cu-K). 



T = 0.42 K [M]. We cannot evaluate the magnetization directly from the HTE 
for temperatures lower than T = 1.5 K due to the lack of convergence. At 
the higher temperatures T = 12.3 K and T = 4.04 K the HTE magnetization 
curves fits are seen to fit the experimental data well. At low temperatures 



60 



M.T. Batchelor, X.-W. Guan, N. Oelkers and Z. Tsuboi 




Figure 18. Magnetization vs magnetic field comparison between theory and experimental 
data 1641 for the compound Cu(Hp)Cl for varying temperature. The TBA with the same coupling 
constants, Jj_ = 14. 5K, Ju = 1.0K and g = 2.1, predicts critical fields at zero temperature which 
coincide with the experimental values at low temperature. The magnetization curves evaluated 
from the HTE give satisfactory agreement with the experimental curves. 



T = 1.6 K and T = 2.44 K, a small discrepancy between the HTE and ex- 
perimental curves is observed. This is mainly because the coupling constants 
J_i_ = 14.5 K and Ju = 1.0 K used in fitting the experimental data of Ref. [59] 
may have some small deviation for fitting the experimental data of Ref. [63] . 
It is observed that the temperature induces a spin-flip in the gapped ground 
state. 



5.4 (C 5 H 12 N ) 2 CuBr 4 



5.4.1 Magnetic suceptibility. The ladder compound (CsH^N^CuBr^ aka 
BPCB [47] . has provided a new ground for understanding microscopic mech- 
anisms in low-dimensional physics. The ladder structure extends along the a- 
axis [47]. The Cu 2+ ions are bridged via the orbital overlap of Br -2 ions along 
the intra- and inter-chain directions. This compound was identified as a spin-^ 
Heisenberg two-leg ladder in the strong coupling regime with J± = 13.3 K and 
J|| = 3.8 K [47]. The gapped phase was observed directly from the experimen- 
tal high field magnetization at low temperature with the typical rounded peak 
susceptibility shape at T = 8K. In fitting the HTE susceptibility to the ex- 
perimental results, we find that the integrable ladder model (|3|) with coupling 
constants J± = 15.4 K, Jy = 1.2 K and g = 2.13 describes this compound well. 
A comparison between the HTE result and experimental susceptibility is given 
in figure fl9l A rounded peak is observed at T 8K. The susceptibility expo- 
nentially decays as the temperature decreases, which indicates a gapped phase. 
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Figure 19. Susceptibility vs temperature comparison between the HTE method and the 
experimental data 1471 for the compound (CsHi2N)2CuBr4. The solid curve is the susceptibility 
evaluated directly from the HTE at H = T. A parameter fit suggests the coupling constants 
Jj_ = 15.4 K and J« = 1.2 K with g = 2.13 and fig = 0.672 K/T. The conversion constant is 

XHTE ~ 0.40615xexp- 

In this case the Troyer et al formula [66] . suggests a gap of 

A ~ 10 K, which coincides with our TBA zero temperature gap A ~ 10.06 K. 

5.4.2 High field magnetization. The high field magnetization has been 
measured at different temperatures [37] ■ The experimental magnetization 
curve at T = 0.7 K reveals that the singlet gapped phase lies in the region 
H < 6.6 T and the saturation magnetization occurs at magnetic fields greater 
than H C 2 = 14.6 T. The TBA magnetization curve is given in the inset of 
figure [20l where the parameter setting is the same as in figure [T9l This curve 
indicates the critical field values H c \ = 7T and H c \ = 14.1 T which are in 
agreement with the experimental result. It is seen that the HTE magnetiza- 
tion curves at T = 1.75 K (solid line) and T = 3.31 K (dashed line) in figure 
[20] fit the experimental curves well. 

5.5 [Cu 2 (C 2 2 )(C 10 H 8 N 2 ) 2 )](N0 3 ) 2 

5.5.1 Magnetic suceptibility. The magnetic properties of the strong cou- 
pling ladder compound [Cu2(C202)(CioHgN2)2)](N03)2 (abbreviated Cu- 
CON) were investigated experimentally in Ref. [161] , The spin-spin exchange 
interactions between Cu 2+ ions along the leg and rung are bridged by NO3 
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Figure 20. Magnetization vs magnetic field comparison between the HTE approach and the 
experimental data |47l for the compound (CsHi2N)2CuBr4. These curves indicate the nature of 
the high magnetic field quantum phase diagram. The same parameters are used as in figure [l9l The 
inset shows the exact TBA magnetization curve obtained at zero temperature. 



and C20^~ , respectively and believed to be antiferromagnetic. The experi- 
mental susceptibility of CuCON has a broad peak at about T = 300 K. Below 
this temperature, the susceptibility smoothly decreases, indicating a large gap. 
Perturbation theory predicts the gap A « Jj_ — Jii th 360 K. The HTE suscepti- 
bility evaluated from the free energy with the coupling constants J± = 520 K, 
J|| = 10 K, g = 2.14 and (i B = 0.672 K/T is shown in figure ED (solid line). An 
overall good fit to the experimental curve is found. The inset of figure I2T1 shows 
the magnetization curves at different temperature. The TBA result predicts 
the critical field values R c \ « 334 T and H C 2 ~ 389 T. These high field values 
appear to be experimentally inaccessible. 



5.6 BIP-BNO 

5.6.1 Magnetic suceptibility. The organic polyradical BIP-BNO has been 
recognised as a strong coupling ladder compound [159. 160] . The intramolecular 
spin-spin exchange interaction within the BIP-BNO molecules is antiferromag- 
netic. Two intermolecular NO groups have the same crystallographic structure 
so that there is a twofold symmetry along the 6-axis. It can be seen to give rise 
to an antiferromagnetic doubled rung interaction. Experimental investigation 
of the magnetic susceptibility reveals a broad peak at T ps 45 K |159[|160j . 
Below this temperature the susceptibility decreases steeply as the tempera- 
ture decreases. Fitting the Troyer et al formula [66], x °£ e~ A / T /VT, below 
T = 25 K gives [159tll60| an energy gap A w 47 K. To fit the susceptibility 
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Figure 21. Susceptibility vs temperature comparison between theory (solid line) and experimental 
data (circles) 11611 for the ladder compound CuCON. The solid line is the susceptibility evaluated 

directly from the HTE at H = T. A parameter fit suggests the coupling constants J± = 520 K 
and Jy = 10 K with g = 2.14 and (/g = 0.672 K/T. The inset shows the HTE magnetization versus 
magnetic field at different temperatures. The conversion constant is Xhte ~ 0.40615xexp- 



and the magnetization, two sets of coupling constants, J± = 73 K, Jm = 17 K 
and J± = 67 K, Jm = 25 K have been suggested [159|I160| . However, the experi- 
mental magnetization curve at temperature T = 1.6 K indicates an energy gap 
A ~ 38.8 T (roughly A ~ 52 K). From perturbation theory, they caculated the 
gap A = Jj_ — Jy = 56 K or 42 K corresponding to the two sets of coupling pa- 
rameters. There is obviously some inconsistency among these fittings. We find 
the parameters Jj_ = 84 K and Jm = 7.0 K with g = 2.0 for the integrable lad- 
der model © in fitting the susceptibility, as shown in figure [22j Here the value 
for the g-factor is fixed in the experiment [159U160] . The solid HTE line shows 
an excellent fit with experimental susceptibility curve. Here the TBA gives 
the gap A ~ 41 T (roughly A ~ 55 K), which is in agreement with the experi- 
mental observation of A ~ 38.8 T. We conclude that their coupling constants 
J± = 73 K and Jm = 17 K seem to be more reasonable for the Heisenberg lad- 
der. In addition, we mention that the conversion constant for this compound is 
twice as large as the one for the other compounds considered here, but of the 
same order as the one for the spin-1 chain compounds [154j . This is mainly 
because of the doubled intermolecular structure along the rungs, effectively 
forming two overlapping spin ladders. 
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Figure 22. Susceptibility vs temperature comparison between HTE theory and 
experiment 159 160 for the strong coupling ladder compound BIP-BNO. The solid line is the 
susceptibility evaluated directly from the HTE at H = 0T. A parameter fit suggests the coupling 
constants = 83 K and Jy = 7.0 K with g = 2.0 and fig = 0.672 K/T. The conversion constant is 
Xhte ~ 2 X 0.406f5xEXP due to the double chain structure. 

6 Thermodynamics of the spin-(i,l) ladder model 

In this section we turn to the calculation of the thermodynamic properties of 
the integrable mixed spin-(i, 1) model (fl"3|) . Various theoretical studies suggest 
that fractional magnetization plateaux exist in mixed spin- (i,l) chains [H5UBH 
EHKHSlETlEHllEniEO] and the mixed spin ladders [82l[8ll[80l[Eni[l9l] . Such frac- 
tional magnetization plateaux have been found in a number of low-dimensional 
magnetic systems, including Shastry-Sutherland systems [76 1 I7 H [7511192] and 
spin ladders [77l[78l[7m [TM[lM[195l . 

6.1 TBA analysis and ground state properties 

A one-third magnetization plateau has been found to exist in the integrable 
mixed spin-(^, 1) ladder for strong rung coupling via the TBA approach |150j . 
Gapped or gapless states appear in turn as the external magnetic field in- 
creases, as can be seen in figure [23] where the whole magnetization curve is eval- 
uated numerically by solving the TBA equations in the different phases pL50j . 
This approach predicts that a two-component massless quantum magnetic 
phase lies in the regime H < H&, with a ferrimagnetic phase appearing for 
H>H c i where the component ^ (recall (|17p ) becomes a physical ferrimag- 

2 

netic ground state. This state forms a one third magnetization plateau which 
extends until the critical field value H C 2- As the magnetic field increases be- 
yond H C 2, the state tjjz becomes involved in the ground state and the ground 
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state becomes a mixture of doublet and quadruplet states. The two compo- 
nents ib[ ^ and ips form an effective XXZ model with massless excitations. 

Other components of the multiplets are gapful by virtue of both the rung cou- 
pling and magnetic field. Eventually, if the magnetic field is greater than H c ^, 
the ground state is the fully polarized state. The critical fields are found from 
the TBA to be Q30] 



H, 



4Jn 



cl 



[gs + - gt)]vB 



3J i — 8J\\ 

H * = rf , /i, " ,i ( 134 ) 
[{9t+9s) - zigs - gt)\t*B 

3J ± + 8J|| 



Hc3 — j, ; : J 



[(gt + g s ) - sigs - gt)]nB 

The magnetization in the vicinity of the critical fields exhibits square root field 
dependent behaviour [150] . 

On the other hand, the first magnetization plateau depends mainly on the 
rung coupling. If | Jn In 2 < J± < J^ F ', where 



J+F = 8J n [\gt + lg s + \{g s - gt)} ^ 135 ^ 
3 [g B + \{g s - gt)] 

this plateau disappears and the critical points H c \ and become close to 
each other. Figure [231 shows the numerical evaluation of the full magnetization 
curve for weak rung coupling J± = 1.3 K. We clearly see that the fractional 
magnetization plateau is closed. This is because the rung interaction can not 
quench the quadruplet. The numerical value is H c \ ~ 0.605 T. This implies 
that for H < H c \ the ground state is the doublet spin liquid phase. The 
quadruplet component tbs is involved in the groundstate for H > H&. There- 

2 

fore the critical point H c \ indicates a quantum phase transition from a two- 
state phase to a three-state phase. Hence for H > H c ±, two Fermi seas, and 
e^ 2 ) , are filled. However, the probability of the component tjj3 quickly increases 

2 

while the probability of the component ip~! quickly decreases as the magnetic 

2 

field increases. It follows that the magnetization increases rapidly in the region 
H c i < H <H C 2, where the critical point H C 2 indicates a phase transition from 
the three-state into the two-state phase. Thus from the theory we predict the 
critical field to be H c2 = 6J±/[(g t +g s ) - \{g s - g t )]^B - H cl « 0.715 T. This 
again is in agreement with the numerical value of = 0.713 T. In the region 
H C 2 < H < H C 3 the two components ipT and ips compete to be in the ground- 
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Figure 23. Magnetization versus magnetic field H for the integrable model l|13| l in the strong 
antiferromagnetic rung coupling regime. The magnetic moment is normalised by the saturation 

value M a = \{g s + gt)fJ-B- The full curve is obtained from the TBA analysis. The coupling 
constants used are |150j J± = 6.0 K, Jn = 0.4 K with Lande g-factors g s = 2.22, gt = 2.09 and 
fig = 0.672 K/T. We predict from the TBA analysis that the one-third saturation magnetization 

plateau opens with the critical fields H cl ft 1.052 T, H c2 ft 5.162 T and H c3 ft 7.394 T which 
coincide with the numerically estimated values. The dashed curve is evaluated from the HTE free 
energy H138I I. Both the TBA and HTE approaches are seen to yield consistent results. In addition, 
the inflection point at H = Hip ft 6.278 T and M ft! 1 is observed from the curves. Hip indicates a 



point of equal probability for the states ij}3 and if>\ 



(-) 



state. If the magnetic field is strong enough, so that H > H c % where H c % is 
given in (|134j) . the reference state ipz becomes the true physical groundstate. 



6.2 HTE approach 

In order to investigate the thermodynamic properties, we again apply the HTE 
scheme which has been set up in section 2J Explicitly, the eigenvalue of the 
QTM (up to constants in the chemical potentials) is given by [131j 

Qi{v - ii) 



+e~^ 2 (j)-(v)cj) + (v) 



Qi(v - %i)Q2(v) 
Qi(v - ±i)Q 2 (v - i) 
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Figure 24. Magnetization versus magnetic field H in the weak antiferromagnetic rung coupling 
regime 150 . The numerical values are the same as for figure l23l but now using the smaller rung 
coupling constant = l.SK < J^ F ~ 2.07K. In this case the fractional magnetization plateau 
vanishes. The TBA critical fields coincide again with the numerically estimated values in the figure. 
The inset shows an enlargement of the magnetization between H c i and H c 2- 
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(136) 



In the above equations, the chemical potentials, m, i = 1, . . . , 6, are chosen as 



Ml 



M2 



f J± - \g s ^Bh + -^J±yi + \{g s ti - gth 



1)2 

2' 
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A»3 
^5 



J± + losfJ-Bh - ^J±y 1 + kidsh' - g t ti + \ f 



J± - \g s ^Bh - -^J±\/ 1 + \(g s h - g t h 



1)2 



Me = {\gt + 9s)vbH, 



(137) 



with h = fisH/J±. These complicated chemical potentials make the HTE 
expansion of the free energy of the integrable model (|13p more involved. As 
a matter of fact, for the applications we have in mind, the third-order term 
in the free energy is rather small. The higher order coefficents do not make a 
significant contribution to the free energy at high temeprature. Therefore, for 
simplicity, we only consider the free energy up to third order, i.e., 



1 



f(T,H) = lnQ\' +c 



C 1 ) J_ „W 
1,0 



+ c 



(1) 

2.0 



+ c 



(1) 

3,0 



J 



+ ... 



(138) 



with the coefficients c[q, c^q and c^q, along with the characters {Q^} for 
su(6), given in Appendix B. 

The exact free energy (|138p allows easy access to the thermodynamic prop- 
erties of the model (|13p via the standard thermodynamic relations given in 
(|129p . Indeed, this approach makes it very easy to predict novel critical be- 
haviour resulting from varying the strength of the rung and leg interactions. 
In figure \23\ the magnetization curve (dashed line) at temperature T = 0.6 K, 
obtained from (|138p , indicates that the ground state remains in a gapless phase 
until the magnetic field exceeds the critical value H c \ ~ 1.5 T. The one-third 
magnetization plateau is also observed under the same parameter setting as 
that used for the TBA analysis. This magnetization plateau vanishes at the 
critical value H C 2 ~ 4.8 T. An ordered ferrimagnetic phase is formed in the 
magnetization plateau region. As the field increases beyond the critical point 
4.8 T, the magnetization increases almost linearly until the whole system is 
fully polarized at the critical value H C 2 ~ 7.8 T. This phase diagram is con- 
sistent with the zero temperature TBA result, also shown in figure [23j This 
behaviour is reminiscent of the integrable spin-| chain |154j . 



The susceptibility curve for this model is given in figure 1251 with the same 
parameters as in figure [231 A gapless antiferromagnetic phase is observed from 
the susceptibility curve. However, the ordered antiferrimagnetic phase results 
in a cusp-like behaviour in the susceptibility. If the magnetization plateau is 
wider, the cusp-like effect becomes more obvious. The inset shows the sus- 
ceptibility for the model with a weaker rung interaction, where the same pa- 
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Figure 25. Susceptibility vs temperature for the model (|13fl with the same parameters as in figure 
1231 The susceptibility curve indicates the antiferromagnetic behaviour in the absence of an energy 
gap. But in this case, the mid-plateau causes a cusp-like feature in the susceptibility curve at 
T fa 2.5 K, similar to what is observed for the integrable spin-^ chain [1 54] . The inset shows the 
susceptibility curve with a weaker rung interaction = 1.3 K, see figure [24] The mid-plateau is 
closed and no cusp-like behaviour is observed in the susceptibility curve. 

rameter settings as figure [Ml are used. As predicted from the TBA analysis, 
the magnetization plateau is closed. Thus the susceptibility shows standard 
antiferromagnetic behaviour without the cusp-like behaviour. 

6.3 Organic ferrimagnetic mixed spin-(^,l) ladder 

To the best of our knowledge, candidates for mixed spin-(i,l) compounds 
appear to be very rare due to their more complicated structure. Nevertheless, 
some organic ferrimagnets may have such mixed structures. One compound 
that we are aware of is the organic ferrimagnet PNNBNO [77 11781179] . which 
has been recognised as a ladder compound with alternating spin-^ and spin- 
1 units, i.e., as two coupled alternating mixed spin chains. PNNBNO forms 
a ferrimagnetic ladder with antiferromagnetic rung and leg interaction. The 
structure along each leg consists of an alternating spin-^ and spin-1 chain, 
which is different from our theoretical model fjl3|> . which has uniform spin-^ 
on one leg and uniform spin-1 on the other. However, at high temperature, the 
thermodynamic properties are entirely dominated by the strong antiferromag- 
netic rung interaction. The integrable model (|13p should therefore be expected 
to adequately describe the compound PNNBNO at high temperatures. In fig- 
ure [26l we give the susceptibility curve evaluated from the exact HTE free 
energy (|138j) . Figure 1261 demonstrates that at high temperatures T > 10 K, 
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Figure 26. Susceptibility vs temperature for the mixed spin-(^, 1) ladder PNNBNO. Circles 
represent the experimental data 1771 . The solid line denotes the susceptibility x ' T evaluated 
directly from the HTE free energy Jl38l with J ± = 220 K and Jy = 0.4 K with \i B = 0.672 K/T 
and g s = gt = 2.03. The conversion constant is Xhte ~ 2 x 0.40615 Xexp due to the double-chain 
structure. The inset shows the magnetization curve predicted from HTE. A large mid-plateau 
forms a ferrimagnetic trimerization phase. 



the theoretical HTE prediction fits the experimental results well. We see that 
the value xT derived from the HTE remains constant in the temperature range 
of 10 — 50 K, while it increases for temperatures over 50 K. This constant range 
reveals the mid-plateau effect - the ground state of the ladder lies in a ferrimag- 
netic state state with S z = |, see the inset of figure EH When T > 50 K, the 
curves indicate antiferromagnetic behaviour. For low temperatures, the HTE 
XT values decrease as the temperature tends to zero. This indicates antiferro- 
magnetic behaviour for the integrable mixed spin-(i, 1) ladder. However, the 
experimental data for the susceptiblity suggests ferromagnetic behaviour at 
low temperature due to the ferrimagnetic structure along the legs. 

In addition, we predict the high temperature entropy at zero field in figure 
[27l It remains constant, i.e., at S = 5.7569 K/J-mol in the temperature range 
10 — 50 K. The entropy increases as the quadruplet enters into the ground 
state from the ferrimagnetic phase (trimerization phase) for T > 50 K. In the 
ferrimagnetic phase, only the doublet component S z = \ occupies the ground 
state. The theoretical prediction for the entropy in this phase is NKb In 2 ~ 
5.76 |77j . which is in excellent agreement with our HTE value S = 5.7569 K/J- 
mol. Below T « 10 K the entropy decreases as the temperature tends to zero. 
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Figure 27. Entropy vs temperature for the the mixed spin-(^, 1) ladder PNNBNO at high 
temperature. The curve evaluated from the HTE with the same constants as in figure [26] for zero 
field. In the temperature region 10 — 50 K, the entropy remains constant at the value S = 5.7569 
K/J-mol, which is consistent with the theoretical value of NKt\n2 a 5.76 K/J-mol |77| . It 
indicates a pure ferrimagnetic trimerization phase with total spin S z = ^. The conversion constant 
is Shte ~ 2 X 4.1535'exp (J/mol-Cu-K), due to the double-chain structure. 



7 Magnetization plateaux for the spin-orbital model 

In this section we discuss extending the analysis to the integrable spin-orbital 
model. 



7.1 Hamiltonians 

The transition metal oxides have been extensively studied in condensed matter 
physics due to their close connections to superconductivity. In particular, the 
electronic and spin properties of the D-electron shells of the transition metal 
ions lead to rich and novel phase transitions. In many metal oxide compounds 
|196lll97lll98pi99ll200| . the orbital order has been observed besides charge 
and spin fluctuations. As a consequence, a great deal of interest |170y20HI202[ 
[M1M[IZI1[IZ21[M1[2QS[I571[I58] in spin-orbital models with both spin and 
orbit double degeneracies arises in this context. In the absence of Hund rule 
coupling, a simple spin-orbital model Hamiltonian is |171|.I172] 




L 

Hs -° = Yl • ^i+i + x ) • T j+ i + y) 

3=1 



(139) 
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where x and y are arbitrary constants, S is a spin-| operator and T is the 
orbital pseudospin-i operator. The model has su{2) (g) su{2) symmetry. 

The phase diagram of this model has been discussed in Refs. |17Hll72|,l205t 
120611207) . In the (x, y)-plane, there are five phases associated with different 
symmetries of the model. In phase I both spin and orbital degrees of freedom 
are in fully polarized ferromagnetic states. In phase II the orbital degrees of 
freedom are in the fully polarized ferromagnetic state while the spin degrees 
of freedom are in the antiferromagnetic ground state and vice versa in phase 
III. Phase IV is a gapped phase. Phase V is a gapless region surrounding the 
origin in the (x,y)-plane. The point (^, ^) on the IV-V phase boundary has 
su(4) symmetry. The simplest su(4) spin-orbital model is given by the per- 
mutation operator. It has a rotation symmetry in the spin S and T spaces. In 
general, this symmetry is broken by Hund's rule coupling or other anisotropic 
hybridizations. In this section, we consider two types of symmetry breaking 
interactions: (1) spin-orbital Ising coupling |157pi58j . and (2) orbital zero-field 
splitting [208] . where the Hamiltonians are respectively given by 

L L 

Hs-o-i = Hs-o-c + Jz ]T S Z T Z - ii bH Y, (9sS z + g t T z ) (140) 

3=1 3=1 
L L 

Hs_o-s = Hs_o-c + A 2 T 3 Z ~ VbH (9sS z + g t T z ) (141) 

i=i i=i 

with 

Hs-o-c = JJ2 ( 2 ^ • 4u + I) [ 2f 3 • Tj+i + I) • (142) 

In the above equations, g s and gt denote Lande factors for spin and orbit. We 
see that for Model 1 (|140p the su(4) symmetry is broken into su(2)®su(2) due 
to the spin-orbital Ising coupling. For Model 2 (]14ip the su{2) symmetry and 
rotation symmetry for the orbit are broken down by the presence of zero-field 
splitting caused by a crystalline field effect. We shall see that these different 
kinds of anisotropy can trigger different phase transitions in the presence of 
an external magnetic field. 

7.2 Ground state properties and fractional magnetization plateaux 

7.2.1 Model 1. In analogy to section[3l the model (|140p can be diagonalised 
by the fundamental basis \S Z ,T Z ), with |1) = |^,^), |2) = i), |3) = 
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| - |, |), |4) = | - i, -i), with energy gSMM 
Ml 1 1 

e = ~ J E rmw-r -^BH( 9s + gt ) Nl 

i=l {Vi ) + 4 

+ - 9t)]N 2 -\\ J z- VBH{g s - g t )]N 3 

+-ti B H(g s + g t )N 4 (143) 

where Aj, i = 1, 2, 3, 4 are the occupation numbers of state \i) and the param- 
eter V{ satisfy the Bethe equations (fTT|) . 



There is no spin singlet state as a result of the different spin and orbit 
Lande factors. The states |2) and |3) have nonzero magnetic moment. The 
Ising coupling J z energetically favours the states |2) and |3). For very strong 
Ising coupling, these two states may form an effective su(2) gapless phase. 
Following the TBA analysis given in section[3l for J z < 4 J In 2, we find that the 
states |1) and |4) are involved in the ground state in the absence of a magnetic 
field. In the presence of a magnetic field, the spin and orbit degrees of freedom 
undergo quite different magnetic ordering due to the Ising coupling. In the 
strong coupling regime J z 3> 4Jln2, the states |1) and |4) are not involved 
in the ground state. When the magnetic field is applied, the spin and orbital 
degrees of freedom are both in an antiferromagnetic state but with opposite 
magnetic moment. Thus the ground state problem is reduced to solving the 
TBA equation 

«<•) W = + Gl - I f_ Q yf^L* (144) 

where the driving term G\ = {g s — Qtj^BH and Q denotes the Fermi boundary. 
The density pW(v) is determined by 

p W( v ) = / 1 v dk. 145) 

P V ^ 2vr v 2 + \ 2rr J_ Q 1 + (v - k) 2 K ' 

If g s > gt, due to Zeeman splitting, the energy level of state |3) becomes 
higher while state |2) becomes lower as the magnetic field increases. As a 
consequence, the magnetization for the spin degrees increases whereas the 
orbital magnetization decreases. If the magnetic field is strong, i.e., H > H&, 
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Figure 28. Magnetization versus magnetic field H in units of saturation magnetization 
Ms = \{9s + S*)Mfl f° r the spin orbital model 11401 . This figure illustrates the magnetization 
plateau obtained from the model via the TBA with J z = 18 K and J = 0.4 K with g s = 2 and 
gt = 1. The inset shows the magnetization for the spin and orbital degrees of freedom. 



with the critical value 

4 T 

Hd = —, r (146) 

VB{g s - gt) 

all sites are occupied by state |2), with no state |3) appearing in the ground 
state. Thus a magnetization plateau with magnetic moment M = ^b{9s — 9t) 
opens up. The spin and orbital degrees of freedom become fully polarized, but 
the orbital moment direction is antiparallel to the magnetic field. Explicitly, 
the total magnetization is M = [Xb9sS z + fiB9tT z , with the spin and orbit 
magnetizations given by 

S z = \- f Q pW(v)dv, T z = ~+ f Q pW(v)dv. (147) 

The spin and orbit magnetization and the total magnetization may be obtained 
by numerically solving (|144p and (|145|) . The results as presented in figure [28l 
On the other hand, the magnetic field can induce the state |1) into the 
ground state. The magnetic plateau vanishes at the critical point 

Hc2 = J -f^~ (148) 
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where the phase transition from a ferrimagnetic phase (with both spin and 
orbital degrees of freedom fully polarized) into an antiferromagnetic phase 
occurs. In this antiferromagnetic phase the spin degrees are fully-polarized 
while the orbital degrees of freedom are in an antiferromagnetic state. We 
notice that in the regime H C 2 < H < the driving term in (1144ft is 

G\ = \j z — fiBdtH. Correspondingly, the magnetization for the spin and 
orbital degrees of freedom is given by S z = \ and T z = — | + J_q (v)dv. In 
the regime ^u^gl — H < ^ci; the driving term is G\ = — \ J z + psgtH, thus 
the magnetization is changed to T z = ^ — JjL p^\v)dv for the orbital degrees 
of freedom. As the magnetic field increases, the state |1) occupies more and 
more sites in the ground state. For strong magnetic fields, i.e., H > H C 3, both 
the spin and orbital degrees of freedom are fully polarized. The critical field 
H c3 is given by 

Hes = 4±^. (149) 

The critical behaviour in the vicinity of the critical points is a Pokrovsky- 
Talapov-type transition as indicated in figure [ 



7.2.2 Model 2. In the above fundamental basis, the nested algebaic Bethe 
Ansatz gives the energy for model 2 in the form 

Mr 1 1 

i=l [v i ) +2 
-[A, + -li B H(g s - g t )]N 2 + {Jt B H{g s - g t )N 3 

-[^-l^Higs-gt)]^. (150) 

Here TVj, i = 1,2,3,4 are the occupation numbers of state \i) and the param- 
eters satisfy the Bethe equations (fTTI) . 

The orbit splitting A z energetically favours the states |2) and |4). For a very 
large orbital splitting, these two states may form an effective su{2) gapless 
phase. Performing the TBA analysis we find that for A z > 2Jln2 the states 
|1) and |3) are gapful in the absence of a magnetic field. Therefore, in this case 
the ground state problem is reduced to the su(2) TBA equation (|144p with the 
driving term G\ = g s fiBH. The density p^(v) is determined by (|145p . In this 
mixed state, the orbital degrees of freedom can be in a fully polarized state, 
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Figure 29. Magnetization versus magnetic field H in units of saturation magnetization 
M s = \ {9s + gt)HB for the spin orbital model d 141 tt . This figure exhibits the magnetization plateau 
obtained through the TBA with A 2 = 8K and J = 0.4 K with g s =2 and gt = 1. The inset shows 
the magnetization for the spin and orbital degrees of freedom. 



while the spin degrees are in an antiferromagnetic phase. Again we assume 
g s > gt, then the energy level of state |4) becomes higher and the energy level 
of state 1 2) becomes lower as the magnetic field increases. As a consequence, 
the magnetization of the spin degrees increases. If the magnetic field is larger 
than the critical value H&, given by 

4 J 

flci = — - (151) 

the state |2) occupies all sites, while the state |4) is not involved in the ground 
state. Here a magnetization plateau with magnetic moment M = ^Bids — gt) 
opens. The spin and orbital degrees of freedom are both fully polarized, as 
seen in the inset of figure I2U1 Explicitly, in the region H < H&, the total mag- 
netization is M = hb9sS z + HB9tT z , with the spin and orbit magnetizations 
given by 

S z = \- f Q pW(v)dv, T z = -\. (152) 

The values of the spin and orbit magnetization and the total magnetization 
may be again obtained by numerically solving (|144p and (|145p . the results of 
which are shown in figure [29l 
For this model the magnetic plateau vanishes for magnetic fields greater 
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than H C 2, due to the involvement of state |1). This critical field has the value 

A Z -4J 

H c2 = . (153) 

VBQt 

At the further critical point H C 2 a phase transition from a ferrimagnetic phase 
into an antiferromagnetic phase occurs. In this phase, the spin state is fully 
polarized, whereas the orbit is in an antiferromagnetic state. As the magnetic 
field increases, state |2) becomes less involved in the ground state. For a strong 
magnetic field, i.e., H > H c ^, both the spin and orbit degrees are fully polar- 
ized. Thus the state |1) becomes the true physical ground state. The critical 
point H C 3 is given by 

ft 2 = ^±^. (154) 

VBQt 

These phase transition points for model 2 (|14ip are clearly indicated in figure 

ES 

We conclude this section by remarking that the thermodynamic properties 
of the spin orbital models are easily accessible from the HTE, along the lines 
discussed earlier in section HI As yet there appears to be a lack of relevant real 
materials, so we do not discuss this further. 



8 Concluding remarks 

In this article we have reviewed a systematic approach for investigating the 
thermal and magnetic properties of the integrable two-leg spin ladder models. 
The thermodynamic Bethe Ansatz (TBA) and High Temperature Expansion 
(HTE) methods used to analyse the ground state properties at zero temper- 
ature and the thermodynamics at finite temperatures have been described in 
detail. The TBA approach outlined in section 3 is best suited to the calcula- 
tion of physical properties at zero temperature. On the other hand, the HTE 
method outlined in section 4 provides an exact expansion for the free energy. 
Although the HTE method involves sophisticated mathematical calculations, 
the resulting expansion for the free energy, e.g., equations (I128P and (|138|) . is 
relatively simple. In practice it suffices to consider only the first few terms in 
the expansion and the relevant thermal and magnetic properties follow almost 
trivially as derivatives to the free energy. In this way the integrable spin-^ 
ladder model defined in section 2.1 can be used to examine the physics of 
a number of strong coupling ladder compounds which have been extensively 
studied in the literature. The two-leg ladder compounds discussed in section 5 
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are (5IAP) 2 CuBr 4 -2H 2 02], Cu 2 (C 5 H 12 N 2 ) 2 Cl4 ISBHES] , (C 5 H 12 N) 2 CuBr 4 
g7], BIP-BNO [T59UTB0] and [Cu 2 (C 2 O 2 )(Ci H 8 N 2 ) 2 )](NO 3 ) 2 OH]. The crit- 
ical magnetic fields derived from the TBA equations are in good agreement 
with the experimental results for the compounds considered. The exact results 
obtained for the susceptibility, the specific heat, the magnetization and the en- 
tropy obtained from the HTE method provide excellent overall agreement with 
the known experimental data. 

In addition to these results, the ground state and thermodynamic properties 
of the integrable mixed spin- (5,1) ladder model defined in section 2.2 have 
been investigated via the TBA and HTE in section 6. The results have been 
used to examine the magnetic properties of the ferrimagnetic mixed spin ladder 
compound PNNBNO [77, 78|[?9"] . Two integrable spin orbital models based on 
an underlying su(4) symmetry were discussed in section 7. It was seen that 
both the single-ion anisotropy and the zero-field orbital splitting can trigger 
a fractional magnetization plateau with respect to different Lande factors for 
the spin and orbital degrees of freedom. The phase diagrams for the two spin 
orbital models were also discussed. 

It is evident from these results that integrable models can provide an alter- 
native way to investigate the thermal and magnetic properties of a number of 
low-dimensional ladder compounds. For arbitrary two-leg spin-(s,5) ladders 
with strong rung coupling, where spins of value s and S reside on different legs, 
the rung interaction terms J± Ylf=i Sj'Sj an d J± Ylj=i(^j'^j) 2 ma y overwhelm 
the spin-spin intrachain interactions and dominate the low temperature be- 
haviour. Therefore it is not surprising that the physics of the spin-(s, S) ladder 
compounds may be captured by the integrable spin-(s, S) ladder models based 
on su(n) symmetry, where n = (2s + 1)(2S + 1). For instance, the integrable 
spin-1 ladder with sophisticated rung interaction may trigger three gapped 
phases with magnetization plateaux at M/M s = 0, ^, 1, which are also found 
to exist in spin chains |209U21(Jll211j . However, the spin-1 ladder compound 
BIP-TENO [2T2l[2l3|[2T4l[2T5] exhibits a magnetization plateau at \ which is 
induced by frustration in the next neighbour interactions along the legs. Ob- 
viously, the integrable spin-1 ladder model does not favor the \ magnetization 
plateau. Nevertheless, it may be used to analyse the fractional magnetization 
plateau for a known spin-1 alternating chain compound [194] . In fact, there 
are relatively few of the more exotic spin ladder compounds. The predictive 
power of this approach has yet to be fully exploited. 

Finally, we emphasize that the gapped phase induced by strong rung in- 
teraction in the spin ladder models is the field-induced singlet ground state 
which is similar to the phase induced by a large single-ion anisotropy for the 
integrable spin chains [153^154] . The transition is of the universal Pokrovsky- 
Talapov type. The universal nature of the critical fields is seen in the full 
H — T phase diagram. For example, the phase diagram for the compounds 
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Figure 30. Phase diagrams for the compounds B5i2aT and Cu(Hp)Cl with spin couplings used in 
sections I5.2l and 15.31 Similar phase diagrams are found using the HTE approach for the other spin 
ladder compounds discussed in section [5] In general the critical point H c 2 approaches the critical 
point H c i as Jii decreases. For Ju = the two critical points merge at zero temperature. 



B5i2aT and Cu(Hp)Cl is shown in Fig. [301 It is anticipated that this phase 
diagram may be further extended into the low temperature region using the 
recently derived nonlinear integral equations for the thermodynamics of the 
su(4) model |216j . In general thermal fluctuations increase the saturation field 
and diminish the energy gap. The slopes of the critical curves indicate that 
the estimated values of H c \ and H C 2 at T = K coincide with the TBA results 
given in table 1. The field-induced gapped phase is significantly different from 
the Haldane valence bond solid state or the dimerization in the spin-Peierls 
transition. The former does not exhibit long range order spin-spin correlations, 
however the latter may. 

The HTE scheme also provides a means of calculating correlation functions 
for integrable models. In particular, there are multiple integral formulae [217] 
(see also, Refs. [2181 122112201 E2H [222]) for the density matrix of the one- 
dimensional Heisenberg model in a magnetic field for any finite temperature. 
The evaluation of the multiple integral formulae is of interest because any 
correlation functions for the model can be written in terms of the density 
matrix. The HTE method has been applied [22311224] for their formulae for 
short ranged case and proved to be useful to evaluate them. In the case of 
zero magnetic field, there is a different way [225 , 226J to calculate the HTE of 
correlation functions. It will be worthwhile to explore correlation functions of 
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the integrable ladder models. 

A particularly exciting development for quantum science in general is that 
fundamental interacting spin systems, traditionally realised in the context of 
condensed-matter based systems like those discussed in this article, may also 
be realised by trapping atoms in optical lattices. This atom-optics based set- 
ting [227] offers a potentially cleaner environment to mimic and study con- 
densed matter physics phenomena. It offers, for example, a way to construct 
a complete toolbox for realising various kinds of spin models [228,229j. It also 
provides a way of realising low-dimensional quantum systems, for which the 
spatial confinement leads to enhanced dynamics and correlations. Integrable 
models, and the methods reviewed in this article for calculating their physi- 
cal properties, will clearly continue to have a bearing on these developments. 
This is because they treat genuine interacting quantum many-body systems for 
which perturbative approaches and mean-field theories often fail. The methods 
discussed here may be generalised to spin systems with open boundary condi- 
tions or impurities [230] . This is particularly important for small-sized systems 
in which quantum fluctuations may be strong enough that boundary effects 
are no longer negligible at low temperature. In general obtaining an exact sys- 
tematic low-temperature expansion for integrable models in analogy with the 
HTE remains an outstanding open problem. So too does the development of 
exact expansion methods for the integrable Bose and Fermi gases. 
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Appendix A: su(4) HTE coefficients 

The first five coefficients in the HTE for the free energy (|128p of the integrable 
spin ladder ([3]) are explicitly given by 



(l) _ 2 Qi 

" Q f 

w _ 3Q< 2) 60j 2 ' 2 , SQf 

" QS" 2 ~ Q?' Q? 3 

(l) _ lOgf 18Qf' 2 80QS 2 ' 3 8QS 3) 24Q< 2) Q< 3) 4Q< 4) 

" 3«3« 2 " Ql 1 ' 4 + 3<3< 1)8 + Qf > 3 ~ <?<» 5 + <?<» 4 

(i) _ 35Qj 2) _ 205QS 2)2 120Qf )3 _ 140Q{ 2)4 55Qj 3) _ 100Q[ 2) Qi 3) 
180Qj 2)2 Qi 3) 45Qj 3)2 15Qj 4) 40Qj 2 >Q{ 4 > 

(i) _ 21Q[ 2) _ 50Qj 2)2 320Q{ 2)3 _ 840Qf )4 4032Qj 2)5 37Qj 3) 
249Qi 2) Qf } 1020Qi L 2)2 Q{ 3) 1344Q( 2)3 Q-; 3) 120Qj 3)2 
378Qj 2) Qj 3)2 34Qj 4) IIVQ^Q^ 336Q? 2 Q? 72Q?Q? 



The Q functions are defined in (|115p . 
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Appendix B: su(6) HTE coefficients 

The first three coefficients in the HTE for the free energy (|138p of the integrable 
mixed spin ladder (fT3|) are explicitly given by 

t-1 n — ^ 



" of 

c 1 =3^-6^ + 3^ 



2 '° " Qi )4 Q? 3 

3,0 3 qS 1)2 q« 4 3 q! 1} qS x) qS 1} 

In turn the characters of the antisymmetric representations of su(6) are ex- 
plicitly written in terms of the chemical potentials (|137p as 

— e -/3(Mi+j"2) _|_ e -/3(/ii+M3) _|_ e -/3(Mi+M4) _|_ g-0(Mi+M5) _|_ e -/3(/ii+Me) 
_|_ e -^(M2+At 3 ) _|_ e -/3(/i 2 +M4) _|_ e -/3(At2+M5) _|_ g-^(M2+Me) _|_ e -/3(ti 3 +LH) 

_|_g-/3(M3+M5) _|_ g-/?(M3+M6) _|_ g-/3(M4+M5) _|_ g-/3(/X4+M6) _|_ g-/3(>5 +A*6 ) 
Q( 3 ) — g-/3(Ml+(U2+M3) _|_ g-/3(Ml+M2+Al4) _|_ e -/3(Mi+M2+Ms) _|_ g-/3(Mi+M2+M(3) 

_|_ e -/3(M2+M3+M4) _|_ g-/3(/i 2 +M3+M5) _|_ g-/3(M2+M3+Me) _|_ e -/3(At2+At4+M."s) 

_|_ e -/3(M2+At4+Me) _|_ g-^(/i 2 +M5+M6) _|_ e -/3(M3+M4+Ats) _|_ e -/3(M3+At4+Me) 

_j_g-/3(M3+M5+M6) _|_ e -/3(M4+/is+Ate) 
q( 4 ) _ g-/3(/ii+M2+M3+/i4) _|_ g-/3(Mi+M2+At3+M5) _|_ e -/3(Ati+M2+M3+Me) 

_|_ e -/3(/Xi+/^ 2 +M4+Ai5) _|_ g-/3(Ml+^2+/i4+M6) _|_ g-/3(/il+/i2+/i5+/i6) 

_|_g-/3(Mi+M3+M4+M5) _|_ g-/3(/-ii+M3+M4+/-ie) _|_ g-/3(Mi+M3+M5+M6) 
_|_g-/ 3 (/- 1 i+/- 1 4+/- 1 5+Me) _|_ g-/3(M2+M3+M4+A*5) _|_ g-/3(M2+M3+M4+Me) 
_)_ e _ / 3 (A 1 2+A 1 3+M5+A 1 6) _|_ g-/3(A 1 2+M4+M5+A 1 6) _|_ g-/3(M3+M4+M5+Me) 
q( 5 ) _ g-/3(Ati+M2+A*3+M4+M5) _|_ g-/3(Mi+M2+M3+At4+Me) _|_ g-^C^+Ms+^+Ms+^e) 

q( 6 ) _ g-/3(Mi+A 1 2+M3+M4+M5+A 1 6) _ q3J±/T 
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Appendix C: Location of Bethe roots and zeros of the largest QTM 
eigenvalue 

In this Appendix the location of Bethe roots and zeros of the largest QTM 
eigenvalue are illustrated for finite systems. 
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Figure CI. (left) Location of roots (a = 1,2,3) of the s«(4) Bethe equations H103I I in the 
complex plane. The parameters are TV = 12, Mi = M2 = M3 = N/2 = 6 and ujy = —0.05, with no 
chemical potentials (fii = for i = 1,2,3,4). Symbols Oi Cj A denote colours 1, 2, 3, respectively. 
Note that the roots of each colour form 6 one-strings, (right) Location of the zeros of the 

corresponding largest QTM eigenvalue T^ a \v) (a = 1,2,3) in the complex plane. The symbols Q, 
□, A denote colours 1, 2, 3, respectively. Note the 'physical strip' around the real axis which is free 

of zeros. 



Integrable quantum spin ladders: comparison between theory and experiment 



89 




0.04 -p. 0.04 0.08 -0.08 -0.04 0.04 0.08 



Re 



Re 



1 1 1 


1 1 1 


A 


A 






1 


1 



Re 



"i 1 i 1 i 1 i 1 r 
o o 

□ □ 

A A 
□ A . . A □ 



O 



A A 

□ □ 



O 



O 



°Re 



o 



Figure C2. (top left) Location of roots v^\a = 1,2,3) of the su(4) Bethe equations 1 11031 in the 
complex plane corresponding to the largest eigenvalue of the QTM for magnetic fields H = 5 T, 9 T 
and 13 T. Symbols o, +, X denote the roots of colour 1, 2, 3, respectively. The parameters are 
N = 12 and T = 1.59 K, with chemical potentials given by equation l|86|l . i.e., fi± = — 13. 3 K, 
fi2 = — 1.4112 H, m = and (14 = —fi2- These values correspond to the compound 
(5IAP)2CuBr4 ■ 2H2O (see figure[4]l with = 13.3 K. The remaining figures are more detailed 

versions of this figure, (top right) Location of the Bethe roots v{ , . Symbols Q, □, A denote 
magnetic fields of H = 5T, 9T and 13 T, respectively, (bottom left) Location of the Bethe roots 



. (bottom right) Location of the Bethe roots v 



(3) 
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Figure C3. (top left) Location of the zeros of the largest QTM eigenvalue T^ a \v) (a = 1,2,3) in 

the complex plane corresponding to the Bethe roots in figure |C21 for the same numerical 
para met ers. This figure is an overlay of the other three figures and corresponds to the top left of 
figure [C2l Symbols o, +, X denote the zeros for colours 1, 2, 3, respectively. Note that the 'physical 
strip' around the real axis is free of zeros for each of the magnetic fields, in accord with Conjecture 
4.1 The remaining panels show the location of zeros of T^ a '(v) for the corresponding Bethe root 

distributions in figure [C2l Symbols 0> I— A denote magnetic fields H = 5T, 9T and 13T, 
respectively, (top right) Colour a = 1. (bottom left) Colour a = 2. (bottom right) Colour a = 3. 



